LOCALISATION AND COMPACTNESS PROPERTIES 
OF THE NAVIER-STOKES GLOBAL REGULARITY 

PROBLEM 



TERENCE TAG 



Abstract. In this paper we establish a number of imphcations 
between various quahtative and quantitative versions of the global 
regularity problem for the Navier-Stokes equations, in the periodic, 
smooth finite energy, smooth H^, Schwartz, or mild categories, 
and with or without a forcing term. In particular, we show that if 
one has global well-posedness in for the periodic Navier-Stokes 
problem with a forcing term, then one can obtain global regular- 
ity both for periodic and for Schwartz initial data (thus yielding a 
positive answer to both official formulations of the problem for the 
Clay Millennium Prize) , and can also obtain global almost smooth 
solutions from smooth data or smooth finite energy data, al- 
though we show in this category that fully smooth solutions are 
not always possible. Our main new tools are localised energy and 
enstrophy estimates to the Navier-Stokes equation that are appli- 
cable for large data or long times, and which may be of independent 
interest. 



1. Introduction 



The purpose of this paper is to estabhsh some imphcations between 
various formulations of the global regularity problem (either with or 
without a forcing term) for the Navier-Stokes system of equations, in- 
cluding the four formulations appearing in the Clay Millennium Prize 
formulation [15] of the problem, and in particular to isolate a single for- 
mulation that implies these four formulations, as well as several other 
natural versions of the problem. In the course of doing so, we also 
establish some new local energy and local enstrophy estimates which 
seem to be of independent interest. 

To describe these various formulations, we must first define properly 
the concept of a solution to the Navier-Stokes problem. We will need 
to study a number of different types of solutions, including periodic 
solutions, finite energy solutions, solutions, and smooth solutions; 
we will also consider a forcing term / in addition to the initial data 
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uq. We begin in the classical regime of smooth solutions. Note that 
even within the category of smooth solutions, there is some choice in 
what decay hypotheses to place on the initial data and solution; for 
instance, one can require that the initial velocity Uq be Schwartz class, 
or merely smooth with finite energy. Intermediate between these two 
will be data which is smooth and in H^. 

More precisely, we define: 

Definition 1.1 (Smooth solutions to the Navier-Stokes system). A 
smooth set of data for the Navier-Stokes system up to time T is a 
triplet {uq, /, T), where < T < cxd is a time, the initial velocity vector 
field Mo : — )■ and the forcing term / : [0, T] x R^ — )■ R^ are 
assumed to be smooth on R^ and [0,T] x R^ respectively (thus, uq 
is infinitely differentiable in space, and / is infinitely differentiable in 
spacetime), and uq is furthermore required to be divergence-free: 

V ■ Mo = 0. (1) 

If / = 0, we say that the data is homogeneous. 

The total energy E{uo, /, T) of a smooth set of data {uq, f, T) is defined 
by the quant it 

E{uo,f,T) := ^(||Mo|Ui(R3) + ||/||LiLi([0,T]xR3))^ (2) 

and (uq,/, T) is said to have finite energy if E{uo,f,T) < oo. We 
define the norm T-L^^uq, f,T) of the data to be the quantity 

V}{uqJ,T) := ||mo||hi(r3) + ||/|U-Hi(R3) < oo 

and say that {uo,f,T) is if 'HMuq, f ,T) < oo; note that the 
regularity is essentially one derivative higher than the energy regularity, 
which is at the level of L^, and instead matches the regularity of the 

initial enstrophy 

u}o{t,x)\'^ dx, 

where uq '■= V x uq is the initial vorticity. We say that a smooth set 
of data (mq, /, T) is Schwartz if, for all integers a,m,k > 0, one has 

sup(l + |x|)'^|V>o(a;)| < oo 
xeR3 

and 

sup (l + |x|)'=|V:9r/(x)|<oo. 
(t,x)e[o,r]xR3 

Thus, for instance, the Schwartz property implies H^, which in turn 
implies finite energy. We also say that [uq, f,T) is periodic with some 

^We will review our notation for spacetime norms such as L^L'^, together with 
sundry other notation, in Section 2 
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period L > if one has uo{x + Lk) = uo{x) and f(t, x + Lk) = f(t, x) 
for all if: G [0,T], x G R'^, and k G Z^. Of course, periodicity is 
incompatible with the Schwartz, H^, or finite energy properties, unless 
the data is zero. To emphasise the periodicity, we will sometimes write 
a periodic set of data (mq, /, T) as {uq, f, T, L). 

A smooth solution to the Navier-Stokes system, or a smooth solution 
for short, is a quintuplet {u,p, Uq, f, T), where (mq, /, T) is a smooth set 
of data, and the velocity vector field u : [0,T] x — )■ and pressure 
field p : [0, T] X R'^ — )■ R are smooth functions on [0, T] x R^ that obey 
the Navier-Stokes equation 

dtu + {u ■ V)m = Am - Vp + / (3) 

and the incompressibility property 

V ■ n = (4) 

on all of [0,T] x R^, and also the initial condition 

u{0, x) = Uo{x) (5) 

for all X G R^. We say that a smooth solution {u,p, Uq, /, T) has finite 
energy if the associated data {uq, f, T) has finite energy, and in addition 
one ha^ 

||'"||Lt°°L2{[0,T]xR3) < OO. (6) 

Similarly, we say that (m, p, Uq, /, T) is if the associated data (uq, /, T) 
is H^, and in addition one has 

Lf^ Hi{[0,T]x-R3) + \\u\\L2Hi{lO,T]xB?) < (7) 

We say instead that a smooth solution {u,p,uq, f,T) is periodic with 
period L > if the associated data (mq, /, T) = (uq, /, T, L) is periodic 
with period L, and if u{t,x + Lk) = u{t,x) for all t G [0,T], a; G R^, 
and k G Z^. (Following [15j, however, we will not initially directly 
require any periodicity properties on the pressure.) As before, we will 
sometimes write a periodic solution {u,p,Uo, f,T) as {u,p,Uo, f,T, L) 
to emphasise the periodicity. 

We will sometimes abuse notation and refer to a solution {u,p, Uq, /, T) 
simply as {u,p) or even u. Similarly, we will sometimes abbreviate a 
set of data (mq,/, T) as {uq, f) or even uq (in the homogeneous case 
/ = 0). 



^Following |15j . we omit the finite energy dissipation condition Vu € 
if L^([0, T] X R'^) that often appears in the Uterature, particularly when discussing 
Leray-Hopf weak solutions. However, it turns out that this condition is actually 
automatic from ([6| and smoothness; see Lemma 8.1 Similarly, from Corollary |ll.l| 
we shall see that the L^H^ condition in ([t]) is in fact redundant. 
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Remark 1.2. In [15], one considerecjj smooth finite energy solutions 
associated to Schwartz data, as well as periodic smooth solutions as- 
sociated to periodic smooth data. In the latter case, one can of course 
normalise the period L to equal 1 by a simple scaling argument. In this 
paper we will be focused on the case when the data (mq, /, T) is large, 
although we will not study the asymptotic regime when T — )> oo. 



We recall the two standard global regularity conjectures for the Navier- 
Stokes equation, using the formulation in [T^ : 

Conjecture 1.3 (Global regularity for homogeneous Schwartz data). 
Let (mo, 0, T) he a homogeneous Schwartz set of data. Then there exists 
a smooth finite energy solution {u,p,uq,0,T) with the indicated data. 

Conjecture 1.4 (Global regularity for homogeneous periodic data). 
Let {uq, 0, T) be a smooth homogeneous periodic set of data. Then there 
exists a smooth periodic solution {u,p,Uq,0,T) with the indicated data. 



In view of these conjectures, one can naturally try to extend them to 
the inhomogeneous case as follows: 

Conjecture 1.5 (Global regularity for Schwartz data). Let {uQ,f,T) 
be a Schwartz set of data. Then there exists a smooth finite energy 
solution {u,p,Uq, f,T) with the indicated data. 

Conjecture 1.6 (Global regularity for periodic data). Let (mq,/, T) 
be a smooth periodic set of data. Then there exists a smooth periodic 
solution {u,p,uo, f,T) with the indicated data. 



As described in [15], a positive answer to either Conjecture 1.3 or Con- 



jecture [L4, or a negative answer to Conjecture [1.5 or Conjecture L6^ 
would qualify for the Clay Millennium Prize. 



However, Conjecture |1.6| is not quite the "right" extension of Con- 
jecture |1.4| to the inhomogeneous setting, and needs to be corrected 
slightly. This is because there is a technical quirk in the inhomoge- 
neous periodic problem as formulated in Conjecture 1.6, due to the 
fact that the pressure p is not required to be periodic. This opens up 
a Galilean invariance in the problem which allows one to homogenise 
away the role of the forcing term. More precisely, we have 



Proposition 1.7 (Elimination of forcing term). Conjecture 1.6 is equiv- 



alent to Conjecture I.4 



The viscosity parameter ly was not normalised in IS] to equal 1, as we are doing 
here, but one can easily reduce to the ^ 1 case by a simple rescaling. 
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We establish this fact in Section [6j We remark that this is the only 
implication we know of that can deduce a global regularity result for the 
inhomogeneous Navier- Stokes problem from a global regularity result 
for the homogeneous Navier-Stokes problem. 



Proposition 1 1 . 7| exploits the technical loophole of non-periodic pressure. 
The same loophole can also be used to easily demonstrate failure of 
uniqueness for the periodic Navier-Stokes problem (although this can 
also be done by the much simpler expedient of noting that one can 
adjust the pressure by an arbitrary constant without affecting ([s])). 
This suggests that in the non-homogeneous case / 7^ 0, one needs an 
additional normalisation to "fix" the periodic Navier-Stokes problem to 
avoid such loopholes. This can be done in a standard way, as follows. 
If one takes the divergence of ([s]) and use the incompressibility Q, one 
sees that that 

Ap = -did,{uiU,) + V ■ f (8) 
where we use the usual summation conventions. If {u,p,Uo, f,T) is a 
smooth periodic solution, then the right-hand side of ([s]) is smooth, 
periodic, and has mean zero. From Fourier analysis, we see that given 
any smooth periodic mean zero function F, there is a unique smooth 
periodic mean zero function A^^F with Laplacian equal to F. We then 
say that the periodic smooth solution {u,p,uo, f,T) has normalised 
pressure if one ha^ 

p = -A-^didjimuj) + A-V ■ /. (9) 

We remark that this normalised pressure condition can also be imposed 
for smooth finite energy solutions (because didj{uiUj) is a second deriv- 
ative of an Ll.{R?) function, and V-/ is the first derivative of an L^(R^) 
function), but it will turn out that normalised pressure is essentially 



automatic in that setting anyway; see Proposition 4.1 



It is well-known that once one imposes the normalised pressure con- 
dition, then the periodic Navier-Stokes problem becomes locally well- 
posed in the smooth category (in particular, smooth solutions are now 
unique, and exist for sufficiently short times from any given smooth 



data); see Theorem 5.1 Related to this, the Galilean invariance trick 



that allows one to artificially homogenise the forcing term / is no longer 



available. We can then pose a "repaired" version of Conjecture 1.6 



Conjecture 1.8 (Global regularity for periodic data with normalised 
pressure). Let {uQ,f,T) be a smooth periodic set of data. Then there 
exists a smooth periodic solution {u,p,uo, f,T) with the indicated data 
and with normalised pressure. 

^Up to the harmless freedom to add a constant to p, this normahsation is equiva- 
lent to requiring that the pressure be periodic with the same period as the solution 
u. 
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It is easy to see that the homogeneous case / = of Conjecture [1.8 
equivalent to Conjecture |1.4[ see e.g. Lemma [4T] below. 



IS 



We now leave the category of classical (smooth) solutions for now, and 
turn instead to the category of periodic mild solutions {u, p, uq, /, T, L). 
By definition, these are functions u, f : [0, T] x R^/LZ^ — )■ R^, p : 
[0, T] X RVLZ^ ^ R, Mo : RVLZ^ R^ with < T, L < cx), obeying 
the regularity hypotheses 

feLrHl{[0,T]x{RyLZ^)) 

u e LTHlr}LlHl{[Q,T] X (rVlz^')) 

with p being given by (|9]), which obey the divergence-free conditions 
(|4]), ([T|, and obey the integral form 

u{t) = e'^uo + f e(*-*')^(-(n ■ V)u -Vp + f){t') dt' (10) 
Jo 

of the Navier-Stokes equation ^ with initial condition ([s]); using the 
Leray projection P onto divergence-free vector fields, we may also ex- 
press (19) equivalently as 

u{t) = e*%o + r e(*-*')^(PE(u, u) + Pf){t') dt' (11) 
Jo 

where B{u, v) is the symmetric bilinear form 

1 

B{u,v)i := --dj{uiVj + UjVi). (12) 

Similarly, we define periodic data to be a quadruplet [uq, f,T, L) 
whose norm 

'H^iuo, f,T, L) := ||mo||//i((r3/lz3)) + ||/IUt°°Hi{(R3/Lz3)) 
is finite, with uq divergence-free. 



Note from Duhamel's formula (20) that every smooth periodic solution 
with normalised pressure is automatically a periodic mild solution. 



As we will recall in Theorem 5.1| below, the Navier-Stokes equation is 
locally well-posed in the periodic category. We can then formulate 
a global well-posedness conjecture in this category: 

Conjecture 1.9 (Global well-posedness in periodic H^)- Let {uq, /, T, L) 
be a periodic set of data. Then there exists a periodic mild so- 
lution {u,p,Uq, f,T, L) with the indicated data. 



We may also phrase a quantitative variant of this conjecture: 
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Conjecture 1.10 {A priori periodic bound). There exists a func- 
tion F : R"*" X R+ X — )■ R+ with the property that whenever 
{u,p, Uq, /, T, L) is a smooth periodic, normalised-pressure solution with 
< T < To < oo and 

H\uoJ,T,L)<A<oc 

then 

||w||Lj»jyi([o,T]xR3/iZ3) < F{A,L,To). 

Remark 1.11. By rescaling, one may set L = 1 in this conjecture with- 
out any loss of generahty; by partitioning the time interval [0, Tq] into 
smaller sub-intervals we may also simultaneously set Tq = 1 if desired. 
Thus, the key point is that the size of the data A is allowed to be 
large (for small A the conjecture follows from the local well-posedness 



theory, see Theorem 5.1). 



As we shall soon see. Conjecture |1.9| and Conjecture |1.10| are actually 
equivalent. 



We now turn to the non-periodic setting. In Conjecture L5, the hy- 
pothesis that the initial data be Schwartz may seem unnecessarily re- 
strictive, given that the incompressible nature of the fluid implies that 
the Schwartz property need not be preserved over time; also, there 
are many interesting examples of initial data that are smooth and fi- 
nite energy (or H^) but not Schwartz. In particular, one can consider 



generalising Conjecture L5 to data that is merely smooth and H^, or 
even smooth and finite energy, rather than Schwartz]^ of Conjecture 



1.5 Unfortunately, the naive generalisation of Conjecture 1.5 (or even 



Conjecture 1.3) fails instantaneously in this case: 

Theorem 1.12 (No smooth solutions from smooth data). There 
exists smooth uo € if^(R^) such that there does not exist any smooth 
finite energy solution {u,p, uo, 0, T) with the indicated data for any T > 
0. 



We prove this proposition in Section 15 At first glance, this propo- 
sition looks close to being a negative answer to either Conjecture L5 
or Conjecture 1.3, but it relies on a technicality; for smooth data. 



the second derivatives of uq need not be square- integrable, and this 
can cause enough oscillation in the pressure to prevent the pressure 
from being (or the velocity field from being Cf) at the initial tim^ 

^We are indebted to Andrea Bertozzi for suggesting these formulations of the 
Navier-Stokes global regularity problem. 

^For most evolutionary PDEs, one can gain unlimited time differentiability at 
t ^ assuming smooth initial data by differentiating the PDE in time (cf. the 
proof of the Cauchy-Kowalesky theorem). However, the problem here is that the 
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t = 0. This theorem should be compared with the classical local ex- 
istence theorem of Heywood [22], which obtains smooth solutions for 
small positive times from smooth data with finite enstrophy, but merely 
obtains continuity at the initial time t = 0. 

The situation is even worse in the inhomogeneous setting; the argument 



in Theorem 1.12 can be used to construct inhomogeneous smooth 
data whose solutions will now be non-smooth in time at all times, not 
just at the initial time t = 0. Because of this, we will not attempt 
to formulate a global regularity problem in the inhomogeneous smooth 
or inhomogeneous smooth finite energy categories. 

In the homogeneous setting, though, we can get around this technical 
obstruction by introducing the notion of an almost smooth finite energy 
solution {u,p,uo, f,T), which is the same concept as a smooth finite 
energy solution, but instead of requiring u,p to be smooth on [0,T] x 
R^, we instead require that u,p are smooth on (0,T] x R^, and for 
each k > 0, the functions Vjiu, dtV^u, V^p exist and are continuous 
on [0,T] X R^. Thus, the only thing that almost smooth solutions 
lack when compared to smooth solutions is a limited amount of time 
differentiability at the starting time t = 0; informally, u is only C^C^ 
at t = 0, and p is only C^C^ at t = 0. This is still enough regularity 
to interpret the Navier-Stokes equation ^ in the classical sense, but 
is not a completely smooth solution. 

The "corrected" conjectures for global regularity in the homogeneous 
smooth and smooth finite energy categories are then 

Conjecture 1.13 (Global almost regularity for homogeneous H^). Let 
(mo,0,T) be a smooth homogeneous set of data. Then there exists 
an almost smooth finite energy solution {u,p, Uq, 0, T) with the indicated 
data. 

Conjecture 1.14 (Global almost regularity for homogeneous finite 
energy data). Let {uq,0,T) be a smooth homogeneous finite energy set 
of data. Then there exists an almost smooth finite energy solution 
{u,p,uq,0,T) with the indicated data. 



We carefully note that these conjectures only concern existence of 
smooth solutions, and not uniqueness; we will comment on some of 
the uniqueness issues later in this paper. 

Another way to repair the global regularity conjectures in these set- 
tings is to abandon smoothness altogether, and work instead with the 

pressure p in the Navier-Stokes equation does not obey an evolutionary PDE, but 
is instead determined in a non-local fashion from the initial data u (see (|9])), which 
prevents one from obtaining much time regularity of the pressure initially. 
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notion of mild solutions. More precisely, define a mild solution 
{u,p, uo, f, T) to be fields u, f : [0, T] x ^ R^, p : [0, T] x R^ ^ R, 
Mo : R^ — )■ R^ with < T < cxD, obeying the regularity hypotheses 

feLrHli[0,T]xR') 

with p being given by (|9|, which obey Q, ([T]), and (10) (and thus 



(11)). Similarly define the concept of data {uo,f,T). 



We then have the following conjectures in the homogeneous setting: 

Conjecture 1.15 (Global well-posedness in homogeneous H^). Let 
(mo,0,T) be a homogeneous set of data. Then there exists a 
mild solution {u,p,Uo,0,T) with the indicated data. 

Conjecture 1.16 {A priori homogeneous bound). There exists 
a function F : R+ x R+ — )• R+ with the property that whenever 
{u,p,Uq,0,T) is a smooth solution with < T < To < oo and 

||Mo|ki(R3) <A<oo 

then 

\\u\\l'^hI{[q,t\x-r?) < F{A,Tq). 
We also phrase a global-in-time variant: 

Conjecture 1.17 {A priori global homogeneous bound). There 
exists a function F : R+ — t- R+ with the property that whenever 
{u,p,uo,0,T) is a smooth solution with 

||'«o|ki(R3) < v4 < OO 

then 

|l'"l|Lt°°i?i([0,T]xR3) < 

In the inhomogeneous setting, we will state two slightly technical con- 
jectures: 

Conjecture 1.18 (Global well-posedness from spatially smooth Schwartz 
data). Let (%,/, T) be data obeying the bounds 

sup(l + |x|)''|V>o(a;)| < oo 

and 

sup (1 + |x|)'=|V^/(x)| < oo 

{t,x)e[0,T]xR3 

for all k,a > 0. Then there exists a mild solution {u,p,uo, f,T) 
with the indicated data. 
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Conjecture 1.19 (Global well-posedness from spatially smooth 
data). Let (mq,/, T) be an set of data, such that 



sup \W^Uo{x) 



< oo 



and 



sup |V"/(a;)| < oo 

{t,x)e[0,T]xK 



for all a > and all compact K . Then there exists a mild solution 
{u,p,Uo, f,T) with the indicated data. 

Needless to say, we do not establish]^ any of these conjectures uncon- 
ditionally in this paper. However, as the main result of this paper, we 
are able to establish the following implications: 

Theorem 1.20 (Implications). 



m 



IV 



fvi 



Conjecture 
Conjecture 



1.9 



1.9 



and Conjecture 1.10 are equivalent, 
implies Conjecture ( and hence also Conjec- 
ture 1.6 and Conjecture 1.4)- 

implies Conjecture 1.19, 
TM 



which is equivalent to 



Conjecture 
Conjecture 



Conjecture 1.19 implies Conjecture 1.13 and Conjecture 1.5 



(and hence also Conjectures 1.3). 



Conjecturell. 13 is equivalent to Conjecture 1.1 4 



Conjecture\1.13, Conjecture 1.15, Conjecture 1.16, and Conjee 



ture \1.1T\ are all equivalent to each other. 



The logical relationship between these conjectures given by the above 
implications (as well as some trivial implications, and the equivalences 
in im ) are displayed in Figure [l] 

Among other things, these results essentially show that in order to 
solve the Navier-Stokes global regularity problem, it suffices to study 
the periodic setting (but with the caveat that one now has to consider 
forcing terms with the regularity of L^H].). 



Theorem 1.20[ i) is a variant of the compactness arguments used in 
(see also (17j. [57]). and is proven in Section [Tj Part (ii) of this theo- 
rem is a standard consequence of the periodic local well-posedness 
theory, which we review in Section |5j In the homogeneous / = case 



^Indeed, the arguments here do not begin to address the main issue in any of 
these conjectures, namely the analysis of fine-scale (and turbulent) behaviour. The 
results in this paper do not prevent singularities from occuring in the Navier-Stokes 
flow; but they can largely localise the impact of such singularities to a bounded 
region of space. 
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Figure 1 . Known implications between the various con- 
jectures described here (i.e. existence of smooth or mild 
solutions, or local or global quantitative bounds in the 
periodic, Schwartz, H^, or finite energy categories, with 
or without normalised pressure, and with or without the 
/ = condition) and also in (the latter conjectures 
and implications occupy the far left column). A positive 
solution to the red problems, or a negative solution to 
the blue problems, qualify for the Clay Millennium prize 
as stated in [15j . 



it is possible to reverse this implication by the compactness arguments 
mentioned previously; see Plj . However, we were unable to obtain this 
converse implication in the inhomogeneous case. Part (iv) is similarly 
a consequence of the non-periodic local well-posedness theory, and 
is also proven in Section |5j 

Part (vi) is also a variant of the results in |31] , with the main new ingre- 
dient being a use of concentration compactness instead of compactness 
in order to deal with the unboundedness of the spatial domain R^, 
using the methods from P, [18], [17] . We establish these results in 
Section m 
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The more novel aspects of this theorem are parts (iii) and (v), which 



we estabhsh in Section [12] and Section [13] respectively. These results 
rely primarily on a new localised enstrophy inequality (Theorem |10.1|_ 
which can be viewed as a weak version of finite speed of propagatioqj 
for the enstrophy ^ |a;(t, x)p dx, where w := V x m is the vorticity. 
We will also obtain a similar localised energy inequality for the energy 
^ /r3 t)ut it will be the enstrophy inequality that is of 

primary importance to us, as the enstrophy is a subcritical quantity 
and can be used to obtain regularity (and local control on enstrophy 
can similarly be used to obtain local regularity). Remarkably, one 
is able to obtain local enstrophy inequahties even though the only a 
priori controlled quantity, namely the energy, is supercritical; the main 
difficulty is a harmonic analysis one, namely to control nonlinear effects 
primarily in terms of the local enstrophy and only secondarily in terms 
of the energy. 

Remark 1.21. As one can see from Figure [TJ the precise relationship 
between all the conjectures discussed here is rather complicated. How- 
ever, if one is willing to ignore the distinction between homogeneous 
and inhomogeneous data, as well as the (rather technical) distinction 
between smooth and almost smooth solutions, then the main implica- 
tions can then be informally summarised as follows: 



• (Homogenisation) Without pressure normalisation, the inhomo- 
geneity in the periodic global regularity conjecture is irrelevant: 
the inhomogeneous regularity conjecture is equivalent to the ho- 
mogeneous one. 

• (Localisation) The global regularity problem in the Schwartz, 
H^, and finite energy categories are "essentially" equivalent to 
each other; 

• (More localisation) The global regularity problem in any of 
the above three categories is "essentially" a consequence of the 
global regularity problem in the periodic category; and 

• (Concentration compactness) Quantitative and qualitative ver- 
sions of the global regularity problem (in a variety of categories) 
are "essentially" equivalent to each other. 



The qualifier "essentially" here though needs to be taken with a grain 
of salt; again, one should consult Figure [T] for an accurate depiction of 
the implications. 



^Actually, in our setting, "finite distance of propagation" would be more ac- 
curate; we obtain an Lj bound for tlie propagation velocity (see Proposition 9.1) 



rather than an bound. 
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The local enstrophy inequality has a number of other consequences, for 
instance allowing one to construct Leray-Hopf weak solutions whose 
(spatial) singularities are compactly supported in space; see Proposi- 
tion [HI 

Remark 1.22. Since the submission of this manuscript, the referee 
pointed out that the partial regularity theory of Cafarelli, Kohn, and 
Nirenberg [4J also allows one to partially reverse the implication in 
Theorem 1.20[iii ), and more specifically to deduce Conjecture 1.8 from 
Conjecture 1.19 We sketch the referee's argument in Remark 12.3 



1.23. Acknowledgements. The author is supported by NSF Research 
Award CCF-0649473, the NSF Waterman Award and a grant from 
the MacArthur Foundation. The author is also indebted to Andrea 
Bertozzi for valuable discussions, and in particular for raising the ques- 
tion of whether the answer to the global regularity problem of Navier- 
Stokes equation is sensitive to the decay hypotheses on the initial ve- 
locity field. We thank Xiaoyi Zhang and the anonymous referee for 
corrections and suggestions. 



2. Notation and basic estimates 



We use A < y, y > X, or X = 0(F) to denote the estimate X < CY 
for an absolute constant C. If we need C to depend on a parameter, 
we shall indicate this by subscripts, thus for instance X <sY denotes 
the estimate X < CgY for some Cg depending on s. We use X ~ F as 
shorthand for X < F < X. 

We will occasionally use the Einstein summation conventions, using 
Roman indices i,j to range over the three spatial dimensions 1,2,3, 
though we will not bother to raise and lower these indices; for instance, 
the components of a vector field u will be Ui. We use di to denote the 
derivative with respect to the i^^ spatial coordinate Xj. Unless otherwise 
specified, the Laplacian A = didi will denote the spatial Laplacian. (In 



Lemma |12.1[ though, we will briefly need to deal with the Laplace- 
Beltrami operator A52 on the sphere S"^.) Similarly, V will refer to 
the spatial gradient V = V^, unless otherwise stated. We use the usual 
notations V/, V ■ m, V x u, for the gradient, divergence, or curl of a 
scalar field / or a vector field u. 

It will be convenient (particularly when dealing with nonlinear er- 
ror terms) to use schematic notation, in which an expression such as 
0{uvw) involving some vector or tensor-valued quantities u,v,w de- 
notes some constant-coefficient combination of products of the compo- 
nents of u, V, w respectively, and similarly for other expressions of this 
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type. Thus, for instance, V x V x m could be written schematically as 
O(V^xi), \u X could be written schematically as 0{uuvv), and so 
forth. 

For any centre Xq G and radius i? > 0, we use B{xq, i?) := {a; G : 
\x — xq\ < R} to denote the (closed) Euclidean ball. Much of our anal- 
ysis will be localised to a ball B{xq, R), an annulus B{xq, R)\B{xo, r), 
or an exterior region Ii^\B{xo, R) (and often xq will be normalised to 
the origin 0). 

We define the absolute value of a tensor in the usual Euclidean sense. 
Thus, for instance, if -u = -Uj is a vector field, then \u\'^ — UiUi, |VMp — 
{diUj){diUj), |V^Mp = {didjUk){didjUk), and so forth. 

If is a set, we use 1e to denote the associated indicator function, 
thus 1e{x) = 1 when x E E and 1e{x) = otherwise. We sometimes 
also use a statement in place of E; thus for instance Ifc^o would equal 
1 if A; 7^ and when A; = 0. 

We use the usual Lebesgue spaces L^{Q) for various domains Q (usu- 
ally subsets of Euclidean space R'^ or a torus R^/LZ^) and various 
exponents 1 < p < oo, which will always be equipped with an obvious 
Lebesgue measure. We often write L*'(Q) as L^^fl) to emphasise the 
spatial nature of the domain O. Given an absolutely integrablc func- 
tion / G L^(R^), we define the Fourier transform / : R^ ^ C by the 
formula 

m I e-^--V(x) dx- 

we then extend this Fourier transform to tempered distributions in the 
usual manner. For a function / which is periodic with period 1, and 
thus representable as a function on the torus R^/Z'^, we define the 
discrete Fourier transform / : — > C by the formula 

/>) := / e-2'^^'=-^/(a;) dx 
7r3/z3 

when / is absolutely integrable on R^/Z'^, and extend this to more 
general distributions on R^/Z^ in the usual fashion. Strictly speaking, 
these two notations are not compatible with each other, but it will 
always be clear in context whether we are using the non-periodic or 
the periodic Fourier transform. 

For any spatial domain (contained in either R^ or R^/LZ^) and any 
natural number A; > 0, we define the classical Sobolev norms ||ii||ijfc(n) 



LOCALISATION AND COMPACTNESS FOR NAVIER-STOKES 



15 



of a smooth function m : ^2 — )■ R by the formula 

k \ V2 

"^'"•'lL2(n) 



Hum ■■= I J] IIv^mI'' 



0=0 

and say that u G H^{^) when ||M||_H'fe(f7) is finite. Note that we do not 
impose any vanishing conditions at the boundary of fl, and to avoid 
technical issues we will not attempt to define these norms for non- 
smooth functions u in the event that Q has a non-trivial boundary. In 
the domain R"^ and for s G R, wc define the Sobolev norm ||ii||ij|(R3) 
of a tempered distribution ti : R^ ^ R by the formula 



1/2 



Strictly speaking, this confiicts slightly with the previous notation when 
A; is a non-negative integer, but the two norms are equivalent up to 
constants (and both norms define a Hilbert space structure), so the 
distinction will not be relevant for our purposes. For s > —3/2, we 
also define the homogeneous Sobolev norm 



1/2 




and let H^(R?), H^(Il^) be the space of tempered distributions with 
finite i7^(R^) or if*(R^) norm respectively. Similarly, on the torus 
R^/Z^ and s G R, we define the Sobolev norm ||ii||ifs(R3/z3) of a dis- 
tribution u : R^/Z^ ^ R by the formula 



F||i/|(R3/Z3) 



again, this conflicts slightly with the classical Sobolev norms if^(R^/Z^), 
but this will not be a serious issue in this paper. We define i7^(R^/Z^) 
to be the space of all distributions u with finite /J|(R^/Z^) norm, and 
i7|(R^/Z'^)o to be the codimension one subspace of functions or distri- 
butions u which are mean zero in the sense that u{0) = 0. 

In a similar vein, given a spatial domain Q and a natural number 
A; > 0, we define C^{fl) to be the space of all k times continuously 
differentiable functions u : D, ^ H whose norm 

k 
j=0 
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is finit^ 

Given any spatial norm associated to a function space de- 
fined on a spatial domain f2, and a time interval J, we can define 
mixed-norms ||M||Lf'Xi(7xQ) functions m : J x f2 — )• R by the formula 

when 1 < p < oo, and 

\\u\\L^x^{ixn) ■■= ess sup \\u{t)\\x,(n), 

assuming in both cases that u{t) lies in X{Q) for almost every Q, and 
then let L^Xx{I x VL) be the space of functions (or, in some cases, distri- 
butions) whose LfXx{I X Q) is finite. Thus, for instance, L'^C'^{I x Q) 
would be the space of functions m : J x f2 — )■ R such that for almost 
every x G /, u{t) : i7 — )■ R is in C^{Q), and the norm 

\\u\\Lf°ci{ixn) ■= ess sup ||u(t)||c2(Q) 
tei 

is finite. 

Similarly, for any natural number > 0, we define C^X^il x Q) to be 
the space of all functions m : / x f2 — )■ R such that the curve t i— )■ u{t) 
from / to Xx{Q) is k times continuously differentiable, and that the 
norm 

k 

\W\\cfX^{IxQ) ■= IIV-'mIIl^oox^c/xq) 

J=0 

is finite. 

Given two normed function spaces X, Y on the same domain (in either 
space or spacetime), we can endow their intersection X (lY with the 
norm 

IklUny := II^^IU + ll^lk- 
For us, the most common example of such hybrid norms will be the 
spaces 

x^i X n) := LTHi^i xn)n lIh^,^\i x n), (13) 

defined whenever J is a time interval, s is a natural number, and Q 
is a spatial domain, or whenever J is a time interval, s is real, and 
n is either R^ or R^Z^ The X' spaces (particularly X^) will play 
a prominent role in the (subcritical) local well-posedness theory for 

^Note that if is non-compact, then it is possible for a smooth function to fail 
to lie in C'^(ri) if it becomes unbounded or excessively oscillatory at infinity. One 
could use a notation such as C^ioc(^) to describe the space of functions that are 
k times continuously differentiable with no bounds on derivatives, but we will not 
need such notation here. 
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the Navier-Stokes equations; see Section [sj The space X° will also 
be naturally associated with energy estimates, and the space with 
enstrophy estimates. 

All of these above function spaces can of course be extended to func- 
tions that are vector or tensor- valued without difficulty (there are mul- 
tiple ways to define the norms in these cases, but all such definitions 
will be equivalent up to constants). 

We use the Fourier transform to define a number of useful multipliers 
on or R^/Z^. On R^, we formally define the inverse Laplacian 
operator A^^ by the formula 

:= ^/(O, (14) 
which is well-defined for any tempered distribution / : — t- R for 



which the right-hand side of (14) is locally integrable. This is for 



instance the case if / lies in k^^ derivative of a function in L;^.(R^) for 
some > 0, or the k^^ derivative of a function in L^.(R^) for some 
k > 1. If / G Ll.(R^), then as is well known one has the Newton 
potential representation 

A"V(x) = ^ [ dy. (15) 

47r Jr. \x - y\ 



Note in particular that (15) implies that if / G LKH^) is supported 
on some closed set K, then A"^/ will be smooth away from K. Also 
observe from Fourier analysis (and decomposition into local and global 
components) that if / is smooth and is either the k^^ derivative of a 
function in L;^(R^) for some > 0, or the k^^ derivative of a function 
in L^(R'^) for some k > 1, then A~^/ will be smooth also. 

We also note that the Newton potential ~ 4^|2_^| is smooth away from 
the diagonal x = y. Because of this, we will often be able to obtain 
large amounts of regularity in space in the "far field" region when \x\ 
is large, for fields such as the velocity field u. However, it will often be 
significantly more challenging to gain significant amounts of regularity 
in time, because the inverse Laplacian A~^ has no smoothing properties 
in the time variable. 

On R^/Z^, we similarly define the inverse Laplacian operator A~^ for 
distributions / : R^/Z^ — )■ R with /(O) = by the formula 
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We define the Leray projection Pu of a (tempered distributional) vector 
field u : — )■ by the formula 

Pu := A"^(V X V X u). 

If u is square-integrable, then Pu is the orthogonal projection of u 
onto the space of square-integrable divergence-free vector fields; from 
Calderon-Zygmund theory we know that the projection P is bounded 
on L^(R^) for every 1 < p < oo, and from Fourier analysis we see that 
P is also ifJ(R^) for every s G R. Note that if u is square-integrable 
and divergence-free, then Pu = u, and we thus have the Biot-Savart 
law 

u = A-\V xu) (17) 

where u := V x u. 

In either R^ or R'^/LZ^, we let e*^ for t > be the usual heat semigroup 
associated to the heat equation ut = An. On R^, this takes the explicit 
form 



for / G L^(R3) for some 1 <p < oo. From Young's inequality we thus 
record the dispersive inequality 

l|e*^/l|L.(R3)<i^"*ll/IU.(R3) (18) 
whenever 1 < p < q < oo and t > 0. 

We recall DuhameVs formula 

u(t) = e(*-*«)^u(to) + /* e(*-*')^(9tu - Au){t') dt' (19) 

whenever u : [to, t] x i7 — )■ R is a smooth tempered distribution, with 
n equal to either R^ or 'R? /7?. 

We record some linear and bilinear estimates involving Duhamel-type 
integrals and the spaces defined in (13), which are useful in the 
local theory for the Navier-Stokes equation: 

Lemma 2.1 (Linear and bilinear estimates). Let [to, ti] he a time inter- 
val, let VL he either R^ or R^/Z^, and suppose that u : [to;^i] x — )• R 
and F : [to,ti] x i7 — )■ R are tempered distrihutions such that 

u{t) = e(*-*")^u(to) + f e(*-*')^F(t') dt' . (20) 

Then we have the standard energy estimat^^ 

||«IU=([to,ti]xn) <s \\u{to)\\Hs(Q) + ||-F||LiH|([to,ii]xn) (21) 

^'^We adopt the convention that an estimate is vacuously true if the right-hand 
side is infinite or undefined. 
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for any s > 0, as well as the variant 

\\u\\x''iito,ti]xn) <s \\u{to)\\Hs(^n] + ll^llL2/fr'([to,ti]xf2) (22) 
for any s > 1. We also note the further variant 

\\u\\x=([to,ti]xn) <s \\u{to)\\H^{n) + \\F\\LiL2(^it^,ti]xn) (23) 
for any s < 3/2. 

We also have the bilinear estimate 

l|v(Mt^)||L4i2([iQ_i^]xf7) < \\u\\xmto,ti]xn)\\v\\xmto,h]x^) (24) 

for any u, v : [to, ti] x — > R, which in particular implies (by a Holder 
in time) that 

\\^{uv)\\L^^L2{[to,U]x-R3) < {ti - to)^^^||'«||xi([to,ii]xR3)||w||xi{[to,ti]xR3). 

(25) 

Proof The estimate^ (g, (g, Q are established in [HI Lemma 



2.1, Proposition 2.2]. The estimate (21) follows from the F = case of 



(21) and Minkowski's inequality. □ 



Finally, we define the Littlewood-Paley projection operators on R^. Let 
(p{^) be a fixed bump function supported in the ball G R^ : |^| < 2} 
and equal to 1 on the ball G R^ : |^| < 1}. Define a dyadic number 
to be a number of the form N = 2^ for some integer k. For each 
dyadic number A^, we define the Fourier multipliers 

iW(0 := ^(e/iV)/(0 
/W(0 := (1 - V{^/N))f{0 

We similarly define P<Ar and P>n- Thus for any tempered distribution 
we have / = J^N^^f in a weakly convergent sense at least, where 
the sum ranges over dyadic numbers. We recall the usual Bernstein 
estimates 

WDTNfWLPCRS) <p,s,D« ^*l|-PAf/llLS(R3), 
II V''PAr/||Lg(R,3) r^k,s ^''||-PAf/||LS(R3), 

3_3„ „ (26) 

II^<Jv/||l«{R3) <p,q Np 1 ||P<7v/||lS(R3), 

3 3 

II^-/v/||lJ(R3) <p,q iV^"^||Pjv/||Lg(R3) 



^"^Strictly speaking, the result in [H] was stated for the torus rather than R'^, but 
the argument works without modification in either domain, after first truncating 
u{to), F to be Schwartz to avoid technicahties at infinity, and usig a standard density 
argument. 
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for all 1 < p < g < oo, s G R, /c > 0, and pseudo-differential operators 
D'^ of order s; see e.g. [l3l Appendix A]. 

We recall the Littlewood-Paley trichotomy: an expression of the form 
PN{{PNifi){PN2f2)) vanishes unless one of the following three scenarios 
holds: 



• (Low-high interaction) N2 ^ Ni ^ N. 

• (High-low interaction) Ni < N2 ^ N. 

• (High-high interaction) N < Ni N2. 



This trichotomy is useful for obtaining estimates on bilinear expres- 
sions, as we shall see in Section [9l 



We have the following frequency- localised variant of (18): 

Lemma 2.2. If N is a dyadic number and f : — ?► R has Fourier 
transform supported on an annulus : |^| ~ A^}, then we have 

l|e*^/||L.(R^) <i^"^exp(-ctAr2)||/||^,(i,3) (27) 
for some absolute constant c > and all 1 < p < q < 00 . 



Proof. By Littlewood-Paley projection it suffices to show that 

\\e'''PNf\\LHm ^ exp(-ctAr2)||/||^,(R3) 

for all test functions /. By rescaling we may set t = 1; in view of 
(18) we may then set A^ > 1. One then verifies from Fourier analysis 
that e*^Pjv is a convolution operator whose kernel has an L^(R^) and 
L;^(R^) norm that are both 0(exp(— cA^^)) for some absolute constant 
c > 0, and the claim follows from Young's inequality. □ 



From the uniform smoothness of the heat kernel we also observe the 
estimate 

\\e^^f\WAK) ^fc,i^,T,p exp(-CTr^)||/||iP(R3) (28) 

whenever 0<t<T, l<p< 00, A; > 0, A^ is a compact subset of R^, 
r > 1, and / is supported on the set {x G R^ : dist(x, A') > r}, and 
some quantity c-r > depending only on T. In practice, this estimate 
will be an effective substitute for finite speed of propagation for the 
heat equation. 



3. Symmetries of the equation 



In this section we review some well-known symmetries of the Navier- 
Stokes fiow, which transform a given smooth solution {u,p,Uo, f,T) 
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to another smooth solution {u,p,uo, f,T), as these symmetries will be 
useful at various points in the paper. 

The simplest symmetry is the spatial translation symmetry 



u{t, x) 
p{t,x) 

f 



u{t, X — Xq) 
p(t,X - Xq) 
Uo{x - Xo) 
f{t,X-Xo) 

T, 



(29) 



valid for any xq G R'^; this transformation clearly maps mild, smooth, 
or almost smooth solutions to solutions of the same type, and also 
preserves conditions such as finite energy, H^, periodicity, pressure 
normalisation, or the Schwartz property. In a similar vein, we have the 
time translation symmetry 



u(t, x) 
p{t,x) 

fit,x) 

f 



u{t + to, x) 
p{t + to,x) 
u(to,x) 
f{t + to,x) 
T-to, 



(30) 



valid for any to ^ [O,^]- Again, this maps mild, smooth, or almost 
smooth solutions to solutions of the same type (and if to > 0, then al- 
most smooth solutions are even upgraded to smooth solutions). If the 
original solution is finite energy or H^, then the transformed solution 
will be finite energy or also. Note however that if it is only the orig- 
inal data that is assumed to be finite energy or H^, as opposed to the 
solution, it is not immediately obvious that the time-translated solu- 
tion remains finite energy or H^, especially in view of the fact that the 
norm (or the enstrophy) is not a conserved quantity of the Navier- 
Stokes flow. (See however Lemma 8.1 and Corollary 11.1 below.) The 
situation is particularly dramatic in the case of Schwartz data; as re- 
marked earlier, time translation can instantly converif^ Schwartz data 
to non-Schwartz data, due to the slow decay of the Newton potential 
appearing in ^ (or of its derivatives, such as the Biot-Savart kernel in 



(17)) 



^^This can be seen for instance by noting that moments such as Jj^3 wi(t, x){x2 ~ 
x\) dx are not conserved in time, but must equal zero whenever u{t) is Schwartz. 
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Next, we record the scaling symmetry 

P{t,x) := ^Pi^,^] 



u,{x) := i^^) (31) 

2 



T := TA , 

valid for any A > 0; it also maps mild, smooth, or almost smooth 
solutions to solutions of the same type, and preserves properties such 
as finite energy, finite enstrophy, pressure normalisation, periodicity, 
or the Schwartz property, though note in the case of periodicity that 
a solution of period L will map to a solution of period XL. We will 
only use scaling symmetry occasionally in this paper, mainly because 
most of the quantities we will be manipulating will be supercritical 
with respect to this symmetry. Nevertheless, this scaling symmetry 
serves a fundamentally important conceptual purpose, by making the 
key distinction between subcritical, critical (or dimensionless) , and su- 
percritical quantities, which can help illuminate many of the results in 
this paper (and was also crucial in allowing the author to discoveij^ 
these results in the first place). 

We record three further symmetries that impact upon the issue of pres- 
sure normalisation. The first is the pressure shifting symmetry 



u(t, x) 
p{t,x) 
tto(x) 
fit,x) 
f 



u(t, x) 

p{t,x) + C{t) 

M^) (32) 

f{t,x) 

T, 



valid for any smooth function C : R — )■ R. This clearly maps smooth, 
or almost smooth solutions to solutions of the same type, and preserves 
properties ssuch as finite energy, , periodicity, and the Schwartz 
property; however, it destroys pressure normalisation (and thus the 
notion of a mild solution). A slightly more sophisticated symmetry in 



^■^The author also found dimensional analysis to be invaluable in checking the 
calculations for errors. One could, if one wished, exploit the scaling symmetry to 
normalise a key parameter (e.g. the energy or a radius parameter r) to equal 
one, which would simplify the numerology slightly, but then one would lose the use 
of dimensional analysis to check for errors, and so we have elected to largely avoid 
the use of scaling normalisations in this paper. 
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the same spirit is the Galilean symmetry 



u{t, x) 


:= u{t, X — v{s] 
Jo 


ds) + v{t) 


p{t,x) 


:=p{t,x— / v{s) 
Jo 

■=uoix) + v{0) 


ds) — X ■ v'{t) 


uo{x) 




f{t,x) 


■= fit,x- / v{s) 
Jo 


ds) 


f 


:=T, 





(33) 



vahd for any smooth function f : R — > R^. One can carefully check 
that this symmetry indeed maps mild, smooth solutions to smooth 
solutions, and preserves periodicity (recall here that in our definition 
of a periodic solution, the pressure was not required to be periodic). 
On the other hand, this symmetry does not preserve finite energy, 
iJ^, or the Schwartz property. It also clearly destroys the pressure 
normalisation property. 

Finally, we observe that one can absorb divergences into the forcing 
term via the forcing symmetry 

■u(t, x) := u{t, x) 

p{t, x) := p{t, x) + q{t, x) 

Uq{x) := Uo{x) (34) 
f{t,x) := f{t,x) + V ■q{t,x), 
f :=T, 



valid for any smooth function P : [0, T] x R^ — R^. If the new forcing 
term / still has finite energy or is still periodic, then the normalisation 



of pressure is preserved. In the periodic setting, we will apply (34) with 



a linear term q{t, x) := x ■ a{t), allowing one to alter / by an arbitrary 



constant a{t). In the finite energy or setting, one can use (34) and 
the Leray projection P to reduce to the divergence-free case V ■ / = 0; 
note though that this projection can destroy the Schwartz nature of 
/. This divergence-free reduction is particularly useful in the case of 
normalised pressure, since ^ then simplifies to 

p = -A~^didj{uiUj). (35) 



One can of course compose these symmetries together to obtain a larger 
(semi)group of symmetries. For instance, by combining (33) and (34) 
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we observe the symmetry 

ft 

v{s) ds) + v{t) 

v{s) ds) 
) 

v{s) ds) + v'{t) 









u(t, x) 


:= u(t, X 


-/'" 






^0 








p{t,x) 


:= p{t, X 


- V 






Jo 


Uo{x) 


■■= Uo{x) 


+ v{0) 


f{t,x) 


■■= fit,x 


-/*" 






^0 


f 


■■=T, 





(36) 



any smooth function f : R — )■ R^. This symmetry is particularly use- 
ful for periodic solutions; note that it preserves both the periodicity 
property and the normalised pressure property. By choosing v{t) ap- 
propriately, we see that we can use this symmetry to normalise periodic 
data (mo, f,T, L) to be mean zero in the sense that 

uo{x) dx = (37) 
and 

f{t,x)dx = (38) 



for all < t < T. By integrating ^ over the torus R^/LZ^, we then 
conclude with this normalisation that u remains mean zero for all times 
< t < T: 

u{t,x)dx = (39) 
The same conclusion also holds for periodic mild solutions. 



4. Pressure normalisation 



The symmetries in (32), (34) can alter the velocity field u and pressure 
p without affecting the data (mq, /, T), thus leading to a breakdown of 
uniqueness for the Navier-Stokes equation. In this section we investi- 
gate this loss of uniqueness, and show that (in the smooth category, 
at least) one can "quotient out" these symmetries by reducing to the 
situation ^ of normalised pressure, at which point uniqueness can be 
recovered (at least in the category). 

More precisely, we show 



Lemma 4.1 (Reduction to normalised pressure). 
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(i) // (m, p, Mo, /, T) is an almost smooth finite energy solution, then 
for almost every time t G [0, T] one has 

p{t,x) = -A-^didj{uiUj){t,x) + ■ f{t,x) + C{t). (40) 

for some bounded measurable function C : [0,T] R. 

(ii) If {u,p,Uq, f,T) is aperiodic smooth solution, then there exists 
smooth functions C : [0,T] — t- R and a : [0,T] — )■ R^ such that 

p(t,x) = -A'^didj{uiUj){t,x) + A'^V-f{t,x)+x-a{t) + C{t). (41) 



In particular, after applying a Galilean transform (33) followed 



by a pressure- shifting transformation (32), one can transform 
{u,p,uo, f,T) into a periodic smooth solution with normalised 
pressure. 

Remark 4.2. Morally, in (i) the function C should be smooth (at least 
for times t > 0), which would then imply that one can apply a pressure- 
shifting transformation ( p2| ) to convert {u,p,Uo, f,T) into a smooth 
solution with normalised pressure. However, there is the technical dif- 
ficulty that in our definition of a finite energy smooth solution, we 
do not a priori have any control time derivatives of u in any L^(R^) 
norms, and as such we do not have time regularity on the component 



A~^didj{uiUj) of (40). In practice, though, this possible irregularity of 
C(t) will not bother us, as we only need to understand the gradient 
Vp of the pressure, rather than the pressure itself, in order to solve the 
Navier-Stokes equations ([3]). 



Proof. We begin with the periodic case, which is particularly easy due 
to Liouville's theorem (which, among other things, implies that the 
only harmonic periodic functions are the constants). We may nor- 
malise the period L to equal 1. Fix an almost smooth periodic solution 
{u,p,uq, f,T). Define the normalised pressure po : [0,T] x R^ — )■ R by 
the formula 

Po := -A-%d,{uiUj) + A- V ■ /. (42) 

As u, f are smooth and periodic, po is smooth also, and from ([s]) one 
has Ap = ApQ. Thus one has 

p = PQ + h 

where h : [0,T] x — )■ R is a smooth function with h{t) harmonic in 
space for each time t. The function h need not be periodic; however, 
from ([3]) we have 

dtu + (m ■ V)m = Au- Vpo -Vh + f. 

Every term aside from Vh is periodic, and so Vh is periodic also. 
Since Vh is also harmonic, it must therefore be constant in space by 
Liouville's theorem. We therefore may write 

h{t,x) = X ■ a{t) + C{t) 
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for some a{t) G and C{t) G R; since h is smooth, a,C are smooth 
also, and the claim follows. 



Now we turn to the finite energy case, thus {u,p,uo, f,T) is now an 
almost smooth finite energy solution. By the time translation sym- 



metry (30) with an arbitrarily small time shift parameter to, we may 



assume without loss of generality that {u,p,Uq, f,T) is smooth (and 



not just almost smooth). We define the normalised pressure pq by (42) 
as before, then for each time t G [0, T] one sees from (fsl) that 



p{t) = Po{t) + h{t) 

for some harmonic function h{t) : R^ — )■ R. As u, f are smooth and 
finite energy, one sees from (42) that po is bounded on compact subsets 
of spacetime; since p is smooth, we conclude that h is bounded on 
compact subsets of spacetime also. From harmonicity, this implies 
that all spatial derivatives V''h are also bounded on compact subsets 
of space time. However, as noted previously, we cannot impose any 
time regularity on po or h because we do not have decay estimates on 
time derivatives of u. 



It is easy to see that h is measurable. To obtain the lemma, it suffices 
to show that h{t) is a constant function of x for almost every time t. 

Let [^1,^2] be any interval in [0,T]. Integrating ^ in time on this 
interval, we see that 



rt2 rt2 

u{t2,x)—u{ti,x)+ {u-'V)u{t, x) dt = / Au{t, x)—\/p{t, x)+f{t, x) dt. 
Jti Jti 

Next, let X : R^ ~^ R be a smooth compactly supported spherically 
symmetric function of total mass 1. We integrate the above formula 
against -^xi^) for some large parameter R, and conclude after some 
integration by parts (which is justified by the compact support of x 
and the smooth (and hence C^) nature of all functions involved) that 



R 3 / u{t2,x)x{^) dx-R ^ u{ti,x)x{^) dx 
Jr3 -K Jr3 rC 



-R- 



t2 



X , 



tl JR3 



u{t, x){u{t, x) ■ Vx)(-^) dxdt = R 



-5 



R' 



t2 



X , 



tl JR3 



u{t, x){Ax){—) dxdt 
R 



+ R-^ 
+ R-^ 



t2 



X 



Vp(t,x)x(— ) dxdt 

tl JR3 R 



t2 



X . 



fit^x)xis) dxdt. 

tl JR3 -K 
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From the finite energy hypothesis and the Cauchy-Schwarz inequahty, 
one easily verifies that 



lim R^^ / u(ti,x)x{-f^) dx = 

?-5-00 Jpj^s R 



lim R'^ / / u{t,x){u{t,x) ■ Vx)(-) dxdt = 
lim R-'^ [ [ u{t,x){Ax){^) dxdt = 
lim R-'^ [ [ f{t,x)x{^) dxdt = 



and thus 



lim i?-^ / / Vp(t,x)x(^) c^xrft = 0. (43) 

iti iR3 



Next, by an integration by parts and (|42j), we can express 

"12 



R-' [ Vpo(t,a;)x(^) dxdt 

Ju JR3 -K 



as 



R-^ [ u^Uj{t, x){VA-^ did jx){^) dxdt 

Jti JRS -f^ 

+ R-^ r ! Mt,x)iVA-'dix){^) dxdt. 



From the finite energy nature of {u, p, uq, f, T) we see that this expres- 
sion goes to zero as i? — i- oo. Subtracting this from (43), we conclude 
that 

Wm R-^ r ! Vh{t,x)x{-) dxdt = 0. (44) 

Jti JR3 R 

'■t2 



The function x i— )■ J^^ Vh{t, x) is weakly harmonic, hence harmonic. 
By the mean-value property of harmonic functions (and our choice of 
x) we thus have 

R-^ [ [ Vh{t,x)x{^) dxdt = [ Vh{t,0) dt 

Jti JR3 Jti 

and thus 

Vh{t,0) dt = 0. 

Since ^1,^2 were arbitrary, we conclude from the Lebesgue differenti- 
ation theorem that Vh{t,0) = for almost every t G [0,T]. Using 
spatial translation invariance (29) to replace the spatial origin by an 
element of a countable dense subset of R^, and using the fact that har- 
monic functions are continuous, we conclude that Vh{t) is identically 
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zero for almost every t G [0, T], and so h{t) is constant for almost every 
t as desired. □ 



We note a useful corollary of Lemma [4.1[ i): 

Corollary 4.3 (Almost smooth solutions are essentially mild). Let 

{u,p,Uo, f,T) be an almost smooth solution. Then {u,p,Uo, f,T) is 
a mild solution, where 



P{t,x) := -A-^didj{uiU,){t,x) + A~'V ■ fit, 



X). 



Furthermore, for almost every t G [0,T], p{t) and p{t) differ by a con- 
stant (and thus Vp = Vp). 



Proof. By Lemma 4.1[ i), Vp is equal to almost everywhere; in 



particular, Vj9 = is a smooth tempered distribution. The claim 
then follows from ([s]) and the Duhamel formula (19). □ 



5. Local well-posedness theory in 



In this section we review the (subcritical) local well-posedness theory 
for both periodic and non-periodic mild solutions. The material 
here is largely standard (and in most cases has been superceded by the 
more powerful critical well-posedness theory); for instance the unique- 
ness theory already follows from the work of Prodi [35] and Serrin |40j . 
the blowup criterion already is present in Leray [31], the local existence 
theory follows from the work of Kato [25], regularity of mild solutions 
follows from the work of Ladyzhenskaya [27], the stability results given 
here follow from the stronger stability results of Chemin and Gallagher 
[6], and the compactness results were already essentially present in 
the author's previous paper [S]. However, for the convenience of the 
reader (and because we want to use the X* function spaces defined 



in (13) as the basis for the theory) we shall present all this theory in 
a self-contained manner. There are now a number of advanced local 
well-posedness results at critical regularity, most notably that of Koch 
and Tataru [26], but we will not need such powerful results here. 

We begin with the periodic theory. By taking advantage of the scaling 



symmetry (31 ) we may set the period L equal to 1. Using the symmetry 



(36) we may also restrict attention to data obeying the mean zero 



conditions (ISTl), (g, thus uq G Hl(R^/Z^)o and / G L^Hl{[0,T] x 
RVZ^)o- 

Theorem 5.1 (Local well-posedness in periodic H^)- Let (mq,/, T, 1) 



be periodic data obeying the mean zero conditions (37), (38). 
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(i) (Strong solution) If {u,p, Uq, /, T, 1) is a periodic mild solu- 
tion, then 

u e C°i/^([0,T] X rVz^). 

In particular, one can unambiguously define u{t) in H].{R? /7?) 
for each t G [0,T]. 

(ii) (Local existence) If 

(||wo||ffi(R3/Z3) + ||/||Lii/i(R3/z3))^T < c (45) 

for a sufficiently small absolute constant c > 0, then there exists 
a periodic mild solution {u,p,Uq, f,T,l) with the indicated 
data with 

lklUi{[0,T]xR3/Z3) < ||'"o||hi(R3/Z3) + ||/||Li/fi{R3/z3) 

and more generally 

I|m||x'=([0,T]xR3/Z3) ~fc,T,||«o!lH|(R3/z3).!l/llLlH|{R3/z3) ^ 

for each k > 1. In particular, one has local existence whenever 
T is sufficiently small depending on 'H}{uq^ f, T, 1). 

(iii) (Uniqueness) There is at most one periodic mild solution 
{u,p,uo, f,T,l) with the indicated data. 

(iv) (Regularity) If {u,p,uq, f,T,l) is a periodic mild solution, 
and (mo,/, T, 1) is smooth, then {u,p,uo, f,T,l) is smooth. 

(v) (Lipschitz stability) Let {u,p,uo, f,T,l) be a periodic mild 
solutions with the bounds < T < Tq and 

||w||xi([0,T]xR3/Z3) < M. 

Let (ug,/',T, 1) be another set of periodic data, and define 
the function 

F{t) := e*^K - «o) + r e(*-*')^(/'(t') - f{t')) dt' . 

Jo 

If the quantity ||-F||xi([o,r]xR3/z3) is sufficiently small depending 
on T, M , then there exists a periodic mild solution {u', p', u^, /', T, 1) 
with 

||m — 'u'||xi([0,T]xR3/Z3) ^T,M ||-^|Ui([0,T]xR3/Z3)- 

Proof. We first prove the strong solution claim (i). The linear solution 

e'^Uo+ re(*-*')^P/(t') dt' 
Jo 



is easily verified to lie in C°i7^([0,T] x R^/Z^), so in view of ([IT]), it 
suffices to show that 

\ii-t')^PB{u{t'),u{t')) dt' 
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also lies in C°if^([0, T] x R'^/Z^). But as u is an mild solution, u lies 
mX\[0,T]xRyZ^), so by (|24|), PB{u,u) lies in L^L2([0, T] x R^/Z^). 
The claim (i) then follows easily from (22). 



|Wo||//i(R3/Z3)- 



Now we establish local existence (ii). Let 5 : 

thus by (45) we have 6^T < c. Using this and (|25|), (|22|) one easily es- 



tablishes that the nonlinear map u defined by 



LiHi(R3/Z3), 



:= e*^Mo + / e(*-*')^(Pfi(M(t'),M(t') + Pf{t')) dt' 
'o 



is a contraction on the ball 

{u G X\[0,T] X RVZ^) : ||m||xi([o,t]xr3/z3) < C6} 

if C is large enough. From the contraction mapping principle we may 
then find a fixed point of $ in this ball, and the claim (ii) follows (the 
estimates for higher k follow from variants of the above argument and 
an induction on k, and are left to the reader). 

Now we establish uniqueness (iii). Suppose for contradiction that we 
have distinct solutions {u,p,uo, f,T,l), {u' ,p' ,uo, f,T,l) for the same 
data. Then we have 

||w||xi([0,T]xR3/Z3), ||M'||xi([0,r]xR3/Z3) < M. 

To show uniqueness, it suffices to do so assuming that T is sufficiently 
small depending on M, as the general case then follows by subdividing 
[0, T] into small enough time intervals and using induction. Subtracting 
(11 ) for u, u' and writing v := u' — u, we see that 

v{t) = f e'^'-'"^^P{2B{u{t'),v{t')) + B{v{t'),v{t')) dt' 



and thus by ( 22 ) 



||'^IUi([o,T]xR3/z3) < MT^/^||t;||xi([o,r]xR3/z3)- 

If T is sufficiently small depending on M, this forces ||f ||xi([o,t]xr3/z3) = 
0, giving uniqueness up to time T; iterating this argument gives the 
claim (iii). 

Now we establish regularity (iv). To abbreviate the notation, all norms 
will be on [0,T] x R^/Z^. As u is an mild solution, it lies in 
X^, hence by (|25|), PB{u,u) lies in L^L^. Applying (|ll|, (|23[, and 
the smoothness of uo,f, we conclude that u G X* for all s < 3/2. 
In particular, by Sobolev embedding we see that u G LfL]^, Vn G 
LiLf n L'^Lf'^ and V^u G L^L^^'^ , hence PB{u,u) G L2ifi([0,T] x 



RVZ^). Returning to (|TT|), (g, we now conclude that u G X^{[0, T] x 
R^/Z^). One can then repeat these arguments iteratively to conclude 
that u G X*^([0,T] x RVZ^) for all k > 1, and thus u G L^C^'([0,T] x 
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RVZ^) for all k>0. From ( [9| w e then have p E L°°C^([0, T] x RVZ^) 
for all A; > 0, and then froni^ we have dtu e L^C*^([0,T] x RVZ^) 
for all k > 0. One can then obtain bounds on dtp and then on higher 
time derivatives of u and t, giving the desired smoothness, and the 
claim (iv) follows. 

Now we establish stability (v). It suffices to establish the claim in 
the short-time case when T is sufficiently small depending only on M 
(more precisely, we take M^T < c for some sufficiently small absolute 
constant c > 0), as the long-time case then follows by subdividing time 
and using induction. The existence of the solution {u',p', u'q, /q, T, 1) is 
then guaranteed by (ii). Evaluating (11) for u,u' and subtracting, and 
setting V := u' — u, we see that 

v{t) = F+ f e(*-*')^P(25(M, v) + B{v, v)){t') dt' 
Jo 



for all t e [0,T]. Applying (p2j, (25), we conclude that 



ll^llxi < \\F\\x^+T'^\\\u\\x^ + \\v\\xi)\\v\\xi 

where all norms are over [to^i] x R-^- Since HwUxi + ll'^llxi is finite, we 
conclude (if T is small enough) that \\v\\x^([o,T]xiiyzS) < ||i^||xi([o,T]xR3/z3), 
as the claim follows. □ 



We may iterate the local well-posedness theory to obtain a dichotomy 
between existence and blowup. Define an incomplete periodic mild 
solution {u,p,uo, f,T~ ,1) from periodic data {uq, f,T^:,l) to be 
fields u : [0, T,) x R^/Z^ -> R^ and v : [0, T*) x R^/Z^ ^ R such that 
for any < T < T^, the restriction {u,p,Uo, f,T, 1) of {u,p,Uo, f,T~, 1) 
to the slab [0,T] x R^/Z'^ is a periodic mild solution. We similarly 
define the notion of an incomplete periodic smooth solution. 

Corollary 5.2 (Maximal Cauchy development). Let {uq, /, T, 1) be pe- 
riodic data. Then at least one of the following two statements hold: 

• There exists a periodic mild solution {u,p,uo, f,T,l) with 
the given data. 

• There exists a blowup time < T^, < T and an incomplete pe- 
riodic mild solution {u,p,Uo, f,T~, 1) up to time T~ , which 
blows up in in the sense that 

lim ||M(t)||iji(R3/z3) = +00. 
We refer to such solutions as maximal Cauchy developments. 

A similar statement holds with "H^ data" and "H^ mild solution" re- 
placed by '^smooth data" and "smooth solution" respectively. 
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Next we establish a compactness property of the periodic flow. 
Proposition 5.3 (Compactness). // (mq""*, /'-"^ T, 1) is a sequence of 



periodic data obeying (37), (38) which is uniformly hounded in 
HliB? /Z^)qxL'^HI{[Q,T]xW /'^ and converges weakl^to (mq, 1), 
and {u,p,uo, f,T,l) is a periodic mild solution with the indicated 
data, then for n sufficiently large, there exists periodic mild solu- 
tions (m*^"\p'^"\ /'^"'^ T, 1) with the indicated data, with u^"^^ con- 
verging weakly in X^{[0,T] x R^/Z^) to u. Furthermore, for any < 
T <T, m'-"-' converges strongly in X^([t,T] x R^/Z^) to u. 

If u"q^ converges strongly in iJ^(R^/Z^)o to uq, then one can set r = 
in the previous claim. 

Proof. This result is essentially in [HI Proposition 2.2], but for the 
convenience of the reader we give a full proof here. 

To begin with, we assume that u*^"-* converges strongly in if^(R^/Z^)o 
to Mo, and relax this to weak convergence later. In view of the stability 



component of Theorem 5.1, it suffices to show that F*^") converges 



strongly in X^([0,T] x R^/Z^) to zero, where 

F(")(t) := e*^(4"^ - Mo) + re(*-*')^P(/("n^') - fit')) dt' . 



We have Mq"'' — Uq converges strongly in if^(R'^/Z^) to zero, while 
_ f converges weakly in L^if^([0,T] x R^/Z^) — > 0, and hence 



strongly in LiLl{\Q,T] x RVZ^). The claim then follows from (22). 



Now we only assume that m*^"-* converges weakly in if^(R^/Z^)o to uq. 
Let < r < T be a sufficiently small time, then from local existence 
(Theorem 5.1|^ii)) we see that m*^"^ and u are bounded in X^([0,r] x 



R /Z ) uniformly in n by some finite quantity M. Writing v 



(n) 



_ then from (|ll|) we have the difference equation 



x;W(t) = pH(t) + / e(*-*')V(5(M,M(")) + 5(M("),M(")))(t') dt'. 



Since Mq"'' — Uq converges weakly in H^R? /7a^) to zero, it converges 



strongly in L^(R^/Z^) to zero too. Using (21) as before we see that 

Fin) 

converges strongly in X°([0,r] x R'^/Z^) to zero. From ( [22| we 
thus have 

ll^^-^llxo < oil) + ||5(m,m("))|L.^-i + ||5(«("\t;("))||^.^-i 



^■^Strictly speaking, we should use "converges in the weak-* sense" or "converges 
in the sense of distributions" here, in order to avoid the pathological (and irrelevant) 
elements of the dual space of L^H^. that can be constructed from the axiom of 
choice. 
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where o(l) goes to zero as n oo, and all spacetime norms are over 
[0, r] X R^/Z^. From the form of B and Holder's inequality we have 

\\u \\l^l^J\v Wl^liW^ IIl2l6 
< Mr 1/^1 It; 1 1x0 
and similarly for B{u,v^"'^), and thus 

ll^^"^IUo<o(l) + Mri/^||t;(")l|xo. 
Thus, for r small enough, one has 

||t;(")|Uo=o(l), 

which among other things gives weak convergence of u^^^ to u in [0, r] x 
R^/Z^. Also, by the pigeonhole principle, one can find times r(") in 
[0, r] such that 

ll^^"H^^"^)lki(R3/Z3)=0(l). 

Using the stability theory, and recalling that r is small, this implies 
that 

||^^("Hr)bi(R3/Z3)=0(l), 

thus u^^\t) converges strongly to u{t). Now we can use our previous 
arguments to extend to all of [0,T] x R'^/Z^ and obtain strong 
convergence in X^{[t,T] x R^/Z^) as desired. □ 



Now we turn to the non-periodic setting. We have the following ana- 



logue of Theorem 5.1 



Theorem 5.4 (Local well-posedness in H^). Let {uq, f,T) be data. 

(i) (Strong solution) If{u,p, uq, /, T, 1) is an mild solution, then 

«eCX'([0,T]xR3). 

(ii) (Local existence and regularity) If 

(lk0b^(R3) + ||/|Lj/^l(R3))^r<C (46) 

for a sufficiently small absolute constant c > 0, then there exists 
a mild solution {u,p,uo, f,T) with the indicated data, with 

lkl|xi([0,T]xR3) < ||mo||/^i(R3) + ||/i|LiHi(R3) 

and more generally 

l|M|Ufe([0,T]xR3) ^fe,ll«0||^.,j,3),||/||^l^fe(j,3),l 

for each k > 1. In particular, one has local existence whenever 
T is sufficiently small depending on H^^uq, f,T) . 

(iii) (Uniqueness) There is at most one mild solution {u,p, Mq, /, T) 
with the indicated data. 
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(iv) (Regularity) If {u,p,uo, f,T,l) is a mild solution, and {uq, f,T) 
is Schwartz, then u andp are smooth; in fact, one has d^u, djp G 
L^H^{[0, T] X R3) for all j,k>0. 

(v) (Lipschitz stability) Let {u,p,Uo, f,T) , {u' ,p' ,Uq, f ,T) be 
mild solutions with the bounds < T < Tq and 

||w||xi([0,T]xR3)7 ||w'||xi([0,r]xR3) < M. 

Define the function 

F{t) := e*^K - «o) + r e(*"*')^(/'(t') - f{t')) dt' . 

Jo 

If the quantity ||-F||l22,2qo.t]xR3) sufficiently small depending 
on T, M , then 

11"" - ^i'||xi([0,T]xR3) ^T,M ||-^||l2l2([0,T]xR3)- 



Proof. This proceeds by repeating the proof of Theorem 5.1 verbatim. 
The one item which perhaps requires some care is the regularity item 
(iv). The arguments from Theorem 5.1 yield the regularity 

u G X''{[0,T] X R3) 

for all A; > without difficulty. In particular, u G L^Hl;{[0,T] x R^) 
for all A; > 0. From ^ and Sobolev embedding one then has p G 
L~i7^([0,T] X R3) for all A; > 0, and then from Q and more Sobolev 
embedding one has dtu G L^if|([0,T] x R'^) for all k > 0. One can 
then obtain bounds on dtp and then on higher time derivatives of u and 
t, giving the desired smoothness, and the claim (iv) follows. (Note that 
these arguments did not require the full power of the hypothesis that 
{uo,f,T) was Schwartz; it would have sufficed to have mq G i/^(R'^) 
and / G C{H^{R^) for all j. A: > 0.) □ 



From the regularity component of the above theorem, we immediately 
conclude that Conjecture 1.19 implies Conjecture L5, which is one half 
of Theorem L20](iv) . 

We will also need a more quantitative version of the regularity state- 
ment in Theorem 15.41 

Lemma 5.5 (Quantitative regularity). Let {u,p,uo, f,T) be an 
mild solution obeying (46) for a sufficiently small absolute constant 
c> 0, and such that 

< M <oo. 



||'«o||ffi(R3) + li; l|LiHfc{R3) 

Then one has 

ll^l|Lt°°Hfc{[T,T]xR3) ^k,T,T,M 1 

for all natural numbers k > 1 and all < t < T . 
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Proof. We allow all implied constants to depend on k, T, M. From 
Theorem 15.11 we have 



lkl|xi([0,T]xR3) < 1 

which already gives the k = 1 case. Now we turn to the k > 2 case. 



From (25) we have 



\\PB{u,u)\\LfLl{[0,T]x-R3) ^ 1 

while from Fourier analysis one has 



JA 



< 1. 



From this and (11), (21) we see that 

||w||x=([r,T]xR3) ^s,T 1 

for all s < 3/2. From Sobolev embedding we conclude 

lkl|L,°°Li2([r,T]xR3) <r 1, 
I|V'u||l2^12([^^T]xR3) 1, 
ll^'"llL^Li^''^([r,T]xR3) ~^ 



and hence 



\PB{u,u)\\L2HU[r,T]x-R3) ^ 1- 



Returning to (11), (23) we now conclude that 

lkllx2([r,r]xR3) 1 

which gives the k = 2 case. One can repeat these arguments iteratively 
to then give the higher k cases. □ 



We extract a particular consequence of the above lemma: 

Proposition 5.6 (Almost regularity). Let {u,p,Uo,0,T) be a homoge- 
neous mild solution obeying (46) for a sufficiently small absolute 



constant c > 0. Then u,p are smooth on [r, T] x for all < t < T; 
in fact, all derivatives ofu,p lie in L'^LI{[t,T] x R-^). If furthermore 
uo is smooth, then {u,p,uo,0,T) is an almost smooth solution. 



Proof. From Lemma 5.5 we see that 

ueLrH',{[T,T]xR') 

for all A; > and < r < T. Arguing as in the proof of Theorem 
5.4 iv) we conclude that u,p are smooth on [r, T] x R^. 
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Now suppose that uq is smooth. Then (since uq is also in iJ^(R^)) 
e*^Mo is smoothp^ on [0,T] x R^, and in particular one has 

r/e*^Mo e L^H^{[0,T] x R^) 

for any smooth, compactly supported cutoff function 7] : R^ — t- R. 
Meanwhile, by arguing as in Lemma 5^ one has 

PB{u, u) G L^LlilO, T] X R^). (47) 



Using (11), (21) one concludes that 

r]u G X'{[0,T] X R3) 

for all cutoff functions rj and all s < 3/2. Continuing the arguments 
from Lemma [5.51 we conclude that 

vPBiu,u)E^Hli[0,T]xR') 



for all cutoffs rj. Using (11), (23) (and using (28), (47), to deal with 
the far field contribution of PB{u,u), and shrinking t] as necessary) 
one then concludes that 

r]u G X\[0,T] X R3) 

for all cutoffs i]. Repeating these arguments iteratively one eventually 
concludes that 

r]u G X^([0,T] X R3) 
for all cutoffs 1], and in particular 

for all > and all compact sets K. By Sobolev embedding, this 
implies that 

u e LrC'^{[0,T] X K) 
for all A; > and all compact sets K. 

We also have m G X^([0, T] x R^), and hence 

ueLrHli[0,T]xR'). 

In particular, 

UiUj e L^Ll{[0,T]xR^) (48) 

and 

UiUj G L'^C^{[0,T] X K) 

for all /c > and compact K. From this and ^ (splitting the inverse 
Laplacian smoothly into local and global components) one has 

peLTCl{\Q,T]xK)- 



^^To obtain smoothness at a point (to, 2:0), one can for instance split Uq into a 
smooth compactly supported component, and a component that vanishes near xq 
but lies in //^^(R-^), and verify that the contribution of each component to e*^Mo is 
smooth at {tQ,xo). 
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inserting this into (|3j) we then see that 

dtue L^C^{[0,T]x K) (49) 
for all A; > and compact K. 

This is a little weaker than what we need for an almost smooth solu- 
tion, because we want V'^m, V'^p, dtV^p to extend continuously down 
to t = 0, and the above estimates merely give LfC^ control on these 
quantities. To upgrade the L"^ control to continuity in time, we first 
observ^^ from (49) and integration in time that we can at least make 
V^u extend continuously to t = 0: 

«eC°C,^([0,T]xif). 

In particular 

UiUj eC^tC^^{[Q,T]x K) (50) 
for all A; > and compact K. 

Now we consider V p in a compact region [0, T] X K. From (g we have 
V^p{t,x) = V^d^dj / TUiU,{t,y) dy. 

Using a smooth cutoff, we spht the Newton potential ^^|^_^| into a "lo- 
cal" portion supported on 5(0, 2i?), and a "global" portion supported 



outside of -8(0, R), where i? is a large radius. From (50) one can verify 
that the contribution of the local portion is continuous on [0, T] x K, 
while from (48) the contribution of the global portion is 0„(l/i?^). 



Sending — oo we conclude that V^p is continuous on [0, T] x K, and 
thus 

pGQ°C,^([0,T]xir) 

for all > and compact K. Inserting this into ([S]) we then conclude 
that 

dtu e C^C^{[0,T] X K) 

for all A; > and compact K, and so we have an almost smooth solution 
as required. □ 

Remark 5.7. Because u has the regularity of L'^H^, we can continue 
iterating the above argument a little more, and eventually get u G 
C^C^{[0,T] X K) and p G ClC^{[0,T] x K) for all A; > and compact 



K. Using the vorticity equation (see (84) below) one can then also get 
u G CfC^{[0,T] X K) as well. But without further decay conditions 
on higher derivatives of u (or of u), one cannot gain infinite regularity 



on u,p,u in time; see Section 15 



^''An alternate argument here would be to approximate the initial data uq by 



Schwartz divergence-free data (using Lemma 12.1) and using a limiting argument 



and the stability and regularity theory in Theorem 5.1 we omit the details. 
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On the other hand, it is possible to use energy methods and the vorticity 
equation (84) to show (working in the homogeneous case / = for 
simphcity) that if Uq is smooth and the initial vorticity Uq := V x uq is 
Schwartz, then the solution in Proposition 5.6 is in fact smooth, with cu 
remaining Schwartz throughout the lifespan of that solution; we omit 
the details. 



As a corollary of the above proposition we see that Conjecture |1.19 



implies Conjecture 1.13 thus completing the proof of Theorem 1.20[iv 



As before, we obtain a dichotomy between existence and blowup. De- 
fine an incomplete mild solution {u,p, uq, /, T~) from data (uq, /, T* 
to be fields u : [0, T*) x -> and v: [0,T^) xR^ ^ R such that for 
any < T < T^, the restriction {u,p,Uo, f,T,l) of {u,p,Uo, f,T~ ,1) 
to the slab [0, T] x R^ is a mild solution. We similarly define the 
notion of an incomplete smooth solution. 

Corollary 5.8 (Maximal Cauchy development). Let {uo,f,T) be 
data. Then at least one of the following two statements hold: 

• There exists a mild solution {u,p,uo, f,T) with the given 
data. 

• There exists a blowup time < T^, < T and an incomplete mild 

solution {u,p,Uq, f,T^) up to time T~ , which blows up in 
the enstrophy norm in the sense that 

lim ||M(t)||Hi(R«) = +00. 



Remark 5.9. In the second conclusion of Corollary |5.8[ more informa- 
tion about the blowup is known. For instance, in the paper [13] it 
was demonstrated that the L^{R^) norm must also blow up (in the 
homogeneous case / = 0, at least). 



6. HOMOGENISATION 



In this section we prove Proposition 1.7 



Fix smooth periodic data {uq, /, T, L); our objective is to find a smooth 
periodic solution {u, p, uq, f, T, L) (without pressure normalisation) with 
this data. By the scaling symmetry pT]) we may normalise the period 
L to equal 1. Using the symmetry (36) we may impose the mean zero 
conditions (37), (38) on this data. 
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By hypothesis, one can find a smooth periodic solution {u, p, uq, 0, T, 1) 



with data (mo,0,T, 1). By Lemma 4.1, and applying a Galilean trans 



form (33 ) if necessary, we may assume the pressure is normalised, which 



in particular makes {u,p,Uq,0,T, 1) a periodic mild solution. 



By the Galilean invariance p3| (with a linearly growing velocity v{t) : = 
2wt), it suffices to find a smooth periodic solution {u,p, Uq, T) (this 
time with pressure normalisation) for the Galilean-shifted data (mq, fw,T), 
where 

and w G is arbitrary. Note that the data (mo,/^,T) continues to 



obey the mean zero conditions (37), (38), and is bounded in if^(R^/Z^)oX 



L^if^([0,T] X R^/Z^)o uniformly in w. We now make a key observa- 
tion: 

Lemma 6.1. If a G R^/Z^ is irrational in the sense that k ■ a ^ 
in R/Z for all k G Z'^\{0}, then fxa converges weakly (or more 
precisely, converges in the sense of spacetime distributions) to zero in 
Lri/i([0,T] xRVZ3)o. 

Proof. It suffices to show that 

/ fxa{t,x)(f){t,x) dxdt — )■ 

for all smooth functions : [0,T] x R'^/Z^ — )■ R. Taking the Fourier 
transform, the left-hand side becomes 

E r e~'-^''''-'^mik)^){-k) dt, 

with the sum being absolutely convergent due to the rapid decrease of 
the Fourier transform of 0(t). Because / has mean zero, we can delete 
the k = term from the sum. This makes k-a non-zero by irrationality, 
and so by the Riemann-Lebesgue lemma, each summand goes to zero 
as A — )■ oo. The claim then follows from the dominated convergence 
theorem. □ 



Let a G R^/Z^ be irrational. By the above lemma, [uq, f\a,T, 1) con- 
verges weakly to (mo,0,T, 1) while being bounded in H^R^ /Z^)q x 
L^if^(R'^/Z^)o. As (mo, 0, T, 1) has a periodic mild solution {u,p, uq, 0, T, 



we conclude from Proposition 5.3 that for A sufficiently large, (mq, /aq, T, 1) 
also has a periodic mild solution, which is necessarily smooth since 
Uq, fxa is smooth. The claim follows. 

Remark 6.2. Suppose that (mq,/, oo,l) is periodic data extend- 
ing over the half- infinite time interval [0, +oo). The above argument 
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shows (assuming Conjecture [1.4 ) that one can, for each < T < oo, 
construct a smooth periodic (but not pressure normahsed) solution 
{u^'^\p^^\uo, f,T,l) up to time T with the above data, by choosing 
a sufficiently rapidly growing linear velocity v^"^"^ = 2w^'^h, applying a 
Galilean transform, and then using the compactness properties of the 
local well-posedness theory. As stated, this argument gives a differ- 
ent solution {nP^\p^^\uQ, /, T, 1) for each time T (note that we do not 
have uniqueness once we abandon pressure normalisation). However, it 
is possible to modify the argument to obtain a single global smooth pe- 
riodic solution (m,p, Mo, /, oo, 1) (which is still not pressure normalised, 
of course), by using the ability in (33) to choose a nonlinear velocity 
f (if:) rather than a linear one. By reworking the above argument, and 
taking v{t) to be a sufficiently rapidly growing function of t, it is then 
possible to obtain a global smooth periodic solution {u^p,Uq, f ^oo^l) 
to the indicated data; we omit the details. 



7. Compactness 



In this section we prove Theorem L20](i), by following the compactness 
arguments of |^4j. By the scaling symmetry (31), we may normalise 
L = 1. 



We ffist assume that Conjecture 1.10 holds, and deduce Conjecture 



L9 
5.2 



Suppose for contradiction that Conjecture |1.9| failed. By Corollary 
there thus exists an incomplete periodic pressure-normalised mild 
solution ('U,p, Mo, /, T~, 1) such that 



lim 



\U 



Hi(R3/Z3) = OO- 



(51) 



By Galilean invariance (36) we may assume that Uq and / (and hence 
u) have mean zero. 

Let (mo"\ f^^\T^,, 1) be a sequence of periodic smooth mean zero data 



converging strongly in if^(R'^/Z^)o 
periodic data {u, f,T^:, 1). For each time < T < T^,, we see from 



X L^Hl{[0,T,] X RVZ^ 



)o to the 



Theorem 5.1 that for n sufficiently large, we may find a smooth solution 







T, 1) with this data, with m'^"-' converging strongly in 



L~i/i([0,T] X R7z3) to u. By Conjecture 1.10, the L'fHl{\{),T] x 



R^/Z^) norm of is bounded uniformly in both T and n, so by 
taking limits as n — )• oo we conclude that ||M(t)||ii-i(R3/z3) is bounded 
uniformly for < t < T.,,, contradicting (51) as desired. 



Conversely, suppose that Conjecture |1.9| held, but Conjecture |1.10 
failed. Carefully negating all the quantifiers, we conclude that there 
exists a time < To < oo and a sequence {u^"'\p''"'\ul^\ f^"'\T^"'\l) 
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of smooth periodic data with < T^") < Tq and /("),T("), 1) 

uniformly bounded in n, such that 

lil^ lkllL?°/fi([0,T(")]xR3/Z3) = oo. (52) 

Using Gahlean transforms (36) we may assume that 4"^/^"^ (and 
hence m*^")) have mean zero. Fro m t he short-time local existence (and 
uniqueness) theory in Theorem 5.1 we see that T^*^^ is bounded uni- 
formly away from zero. Thus by passing to a subsequence we may 
assume that T*^"^ converges to a limit T^, with < T^, < Tq. 

By sequential weak compactness, we may pass to a further subsequence 
and assume that for each < T < T,, (4"\ /("\ T, 1) converges weakly 
(or more precisely, in the sense of distributions) to a periodic limit 
{uo, f,T, 1); gluing these limits together one obtains periodic data 
(mo, /, T^,, 1), which still has mean zero. By Conjecture 1.9, we can then 
find a periodic mild solution {u,p, uq, f, T^,, 1) with this data, which 
then necessarily also has mean zero. 



By Theorem 5J^ and Proposition 5^, we see that for every < r < 
T < r„ converges strongly in L^i7^([r,r] x RVZ^) to u. In 
particular, for any < T < T,,, one has 

limSUp ||m(")(T)||^1(r3/z3) < \\u\\Lf'Hl{[0,T,]x-R:^/Zi) < oo. 

Taking T sufficiently close to T^, and then taking n sufficiently large, 
we conclude from Theorem 15. II that 

limsup||M(")||i^^i([5._r(„)]xR3/z3) < oo; 

also, from the strong convergence in L^if^([r, T] x R'^/Z'^) we have 
limsup ||m(")||lc«^i([^_j-]xr3/z3) < oo 

n—^oo 

for any < r < T, and finally from the local existence (and uniqueness) 
theory in Theorem 5J^ one has 

limsup ||m^"^||looj^i([o^^]xr3/z3) < oo 

n—^ca 

for sufficiently small r. Putting these bounds together, we contradict 
(52), and the claim follows. 

Remark 7.1. It should be clear to the experts that one could have 
replaced the regularity in the above conjectures by other subcrit- 
ical regularities, such as H'' for A; > 1, and obtain a similar result to 
Theorem 1.20[ i). 
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As remarked previously, the homogeneous case / = of Theorem 
1.20[ i) was estabhshed in |33]. We recall the main results of that paper. 
We introduce the following homogeneous periodic conjectures: 

Conjecture 7.2 (^4 priori homogeneous periodic bound). There 
exists a function F : x R"*" x R+ — )■ R+ with the property that when- 
ever {u,p,Uo,0,T, L) is a smooth periodic, homogeneous normalised- 
pressure solution with < T < Tq < oo and 

n\uo,0,T,L) <A<oo 

then 

\\u\\L^HlX[0,T]x-R3/LZ'i) < F{A,L,To). 

Conjecture 7.3 {A priori homogeneous global periodic bound). 
There exists a function F : R"*" x xR+ — R+ with the property that 
whenever {u,p, uq, 0, T, L) is a smooth periodic, homogeneous normalised- 
pressure solution with 

'H\uo,0,T,L) <A<oo 

then 

ll'"||Lj°°Hi([o,r]xR3/iz3) < F{A,L). 

Conjecture 7.4 (Global well-posedness in periodic homogeneous H^). 
Let {uo,0,T, L) be a homogeneous periodic set of data. Then there 
exists a periodic mild solution {u,p,Uq,0,T, L) with the indicated 
data. 

Conjecture 7.5 (Global regularity for homogeneous periodic data 
with normahsed pressure). Let {uo,0,T) be a smooth periodic set of 
data. Then there exists a smooth periodic solution {u,p,Uo,0,T) with 
the indicated data and with normalised pressure. 



In \AA\ Theorem 1.4] it was shown that Conjectures 1.4, 7.2, 7.3 



equivalent. As implicitly observed in that paper also, Conjecture 1.4 



are 
is 



equivalent to Conjecture 7.5 (this can be seen from Lemma 4.1), and 



from the local well-posedness and regularity theory (Theorem 5.1 



or 



[441 Proposition 2.2]) we also see that Conjecture 7.5 is equivalent to 



Conjecture 7.4 



8. Energy lggalisatign 



In this section we establish the energy inequality for the Navier-Stokes 
equation in the smooth finite energy setting. This energy inequality is 
utterly standard (see e.g. |38]) for weaker notions of solutions, so long 
as one has regularity of LfH^, but (somewhat ironically) requires more 
care in the smooth finite energy setting, because we do not assume a 
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priori that smooth finite energy solutions lie in the space L^H^. The 
methods used here are local in nature, and will also provide an energy 



localisation estimate for the Navier-Stokes equation (see Theorem 8.2). 
We begin with the global energy inequality. 

Lemma 8.1 (Global energy inequality). Let {u,p,uo, f,T) be a finite 
energy almost smooth solution. Then 

||M||Lt°°L2{[0,T]xR3) + II VM||i2^2([o,T]xR3) < E{uo, f , Tf'"^. (53) 

In particular, u lies in the space X^{[{],T] x R^). 

Proof. To abbreviate the notation, all spatial norms here will be over 
R3. 



Using the forcing symmetry (34), we may set / to be divergence-free, 
so in particular by Corollary |4.3| we have 

Vp(t) = Vp{t) (54) 

for almost all times t, where 

p= -A'^didj{uiUj). (55) 

As {u,p,Uo, f,T) is finite energy, we have the a priori hypothesis 

||m||lj»l2([0,T]xR3) < A 

for some A < oo, though recall that our final bounds are not allowed to 
depend on this quantity A. Because u is smooth, we see in particular 
from Fatou's lemma that 

\Ht)\\Li<A (56) 

for all t e [0,T]. 

Taking the inner product of the Navier-Stokes equation ^ with u and 
rearranging, we obtain the energy density identity 

dtihu\^) + u ■ V(^|Mn = A(^|Mn - |Vm|2 -u-Vp + u-f. (57) 

We would like to integrate this identity over all of R'^, but we do not yet 
have enough decay in space to achieve this, even with the normalised 
pressure. Instead, we will localise by integrating the identity against a 
cutoff 7]'^, where r]{x) := x{ ^^^~'^ ); and x '■ ^ ^ R-"*" is a fixed smooth 
function that equals on [0, -|-oo] and 1 on [— oo, —1], and < r < R/2 
are parameters to be chosen later. (The exponent 4 is convenient for 
technical reasons, in that r]^ and V{r]^) share a large common factor 
rj^, but it should be ignored on a first reading.) Thus we see that rj^ 
is supported on the ball 5(0, R) and equals 1 on B{0, R — r), with the 
derivative bounds 

V^r/ = 0{r-^) (58) 
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for j = 0, 1, 2. We define the locahsed energy 



Ej^4.{t) := / -\u\'^{t,x)r]^{x) dx. 



(59) 



Clearly we have the initial condition 

E,40)<Eiuo,f). (60) 

Because r/^ is compactly supported and u is almost smooth, _E^4 is C/, 
and we may differentiate under the integral sign and integrate by parts 
without difficulty; using (54), we see for almost every time t that 

dtE^. = -Xi + X2 + Xs + X4 + X5 



where Xi is the dissipation term 



Xi := 

X2 is the heat flux term 

X2 

X3 is the transport term 

X3: 

X4 is the forcing term 



Vm|^?7^ dx = ||?7"'Vm||^2; 



R3 



\ I |«rA(r/^) dx, 



(61) 
(62) 



R3 



\u\'^u ■ rfVr] dx, 



X4 



u ■ frf dx. 



R3 



and X5 is the pressure term 

X5 := 4 / uprfVr] dx. 

JRS 

The dissipation term Xi is non-negative, and will be useful in control- 
ling some of the other terms present here. The heat flux term X2 can 
be bounded using (56) and (58) by 



X, < 



A' 



2 



so we turn now to the transport term X3. Using Holder's inequality 

(63) 



and (58), we may bound 



and thus by (56) and Sobolev embedding 

Xs<—\Muv')\\%^ 



By the Leibniz rule and (62), (56), (58) one has 

A 
r 



|v(V)IIl|<x^/' 
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and thus 

Xs<—x'/' + —. 

Now we move on to the forcing term X4. By Cauchy-Schwarz, we can 
bound this term by 

X4 < Eli'ait) 

where a{t) := ||/(t)||L2(_B(o,i?))- Note from ^ that 

r-r 

a{t) dt<E{uoJ,TY'\ (64) 



Now we turn to the pressure term X5. From (55) we have 



X^= f 0{u{A-^\/\uu))TfVr]). 

JR3 



We will argue as in the estimation of X4, but we will first need to move 
the 7f weight past the singular integral A^^V^. We therefore bound 
= Xs^i + X5,2 where 



and 



X5,i= / 0(m(A-V(uV))V77) 

JR3 



where [A, B] := AB — BA is the commutator, and rj^ is interpreted as 
the multiplication operator yf : u ^ rfu. For Xg^i, we apply Holder's 
inequality and ( 58 ) to obtain 

X5,i<^||«|U.||A-V2(nV)||^^_ 

The singular integral A~^V^ is bounded on L^, so it may be discarded; 
applying Holder's inequality again we conclude that 

V <- ^ W l|3/2|| 2||3/2 

^5.1 ~ -Mil M 



This is the same bound (63) used to bound X3, and so by repeating 
the X3 analysis we conclude that 



As for X5 2, we observe from direct computation of the integral kernel 
that when r = 1, [A~^V^, y^] is a smoothing operator of infinite order 
(cf. [25]), and in particular 

ll[A-V^r/VlU|< ll/IUi 
in the r = 1 case. In the general rescaling argument then gives 

\\[A-'v\m\Li<^\\f\\Li. 
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Applying Holder's inequality and (56) we conclude that 



X, , < 



'5.2 ~ ^5/2' 

Putting all the estimates together, we conclude that 
By Young's inequality we have 



2 

and 



^ A^ 



and so we obtain 

A^ A^ 

and hence for almost every time t 



dtE^. +Xi<^ + ^ + Eli^ait) (65) 



dt{E,. + E{uoJ,T)f" < EiuoJ^T)-'/' i^— + — j +a(t). 

By the fundamental theorem of calculus, (64) and (60), we conclude 
that 

EAt)'^' < E{uoJ.Tf/' + E{uoJ.T)-'" + ^) ^ 

for all t G [0, T] and all sufficiently large i?; sending r, R oo and 
using the monotone convergence theorem we conclude that 

||m|U?°l|{[o,t]xr3) < E{uo,f,TY^^. 
In particular we have 

EAt)<E{uo,f,T) 



for all r, R; inserting this back into (65) and integrating we obtain that 

T / A'^ A^\ 

Xiit)dt<^-^ + ^jT + EiuoJ,T). 

Sending r, i? — )■ oo and using monotone convergence again, we conclude 
that 

l|VM||L2i2([o,T]xR3) < E{uo,f,TY^'^ 

and Lemma |8JJ follows. □ 
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We can bootstrap the proof of Lemma 8A_ as follows. A posteriori, we 
see that we may take A < E{uo, f, T)^/^. If we return to (65), we may 
then obtain 

where e > is an arbitrary parameter which we will optimise later. 
From the fundamental theorem of calculus we then have 



L}Ll, 



where the L^Ll. norm is over [0, T] x i?(0, i?); optimising in e, we con- 
clude that 

Inserting this back into (65) and integrating, we also conclude that 

[x.ii)dt< (b,,(0)./^+ + 



Applying spatial translation invariance (29) to move the origin from 



to an arbitrary point Xq, we conclude an energy localisation result: 

Theorem 8.2 (Local energy estimate). Let {u,p,uo, f,T) be a finite 
energy almost smooth solution with f divergence-free. Then for any 
Xo G and any < r < R/2, one has 

\W\\L^L2{[0,T]xB{xo,R~r)) + II V'u||i2^2([o,T]xB(a:o,iJ-r)) 

\uo\\l2(B(xq,R)) + 



< 



LjLlilO,T]xB{xo,R)) 



(66) 



Remark 8.3. One can verify that the estimate (66) is dimensionally 



consistent. Indeed, if L denotes a length scale, then r, i?, £'(mo, /) have 
the units of L, T has the units of L^, u has the units of and all 



terms in (66) have the scaling of L^^^. Note also that the global energy 



estimate (8.1) can be viewed as the limiting case of (66) when one 



sends r, R to infinity. 



Remark 8.4. A minor modification of the proof of Theorem 8.2 allows 
one to replace the ball B{xo,R) by an annulus 

B{xo,R')\Bixo,R) 

for some < R < R' with r < (R' - R)/2,R/2, with smaller ball 
B{xq, R — r) being replaced by the smaller annulus 

B{xo,R' -r)\B{xo,R + r). 
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The proof is essentially the same, except that the cutoff rj has to be 
adapted to the two indicated annuli rather than to the two indicated 
balls; we omit the details. Sending R' ^ oo using the monotone con- 
vergence theorem, we conclude in particular an external local energy 
estimate of the form 

\W\\L^Ll{[0,T]x{-R3\B{xo,R+r))) + II Vm|| j^2^2Qo_T]x(R3\B(a:o,H+r))) 

^ \\uo\\L'i{R3\B{xo,R)) + \\f\\LlLl{[0,T]x{R^\B{xo,R))) 

EK,/,T)V^TV^ E{u,J\Tfl^T'/^ 

2 

(67) 

whenever < r < R/2. 

Remark 8.5. The hypothesis that / is divergence-free can easily be 



removed using the symmetry (34), but then / needs to be replaced by 



Pf on the right-hand side of (66). 



Remark 8.6. Theorem 8.2 can be extended without difficulty to the 
periodic setting, with the energy E{uq, /, T) being replaced by the pe- 
riodic energy 

El{uoJ,T) := -(||Mo||l2(R3/LZ3) + \\f\\LlLl{[0,T]x-R? /LZ^)f ^ 

as long as the radius R of the ball is significantly smaller than the 
period L of the solution e.g. R < L/100. The reason for this is that 
the analysis used to prove Theorem |8.2| takes place almost entirely 
inside the ball B{xo,R), and so there is almost no distinction between 
the finite energy and the periodic cases. The only place where there is 
any "leakage" outside of B{xo, R) is in the estimation of the term X5 2, 
which involves the non-local commutator [A~^V^, r]^]. However, in the 
regime R < L/100 one easily verifies that the commutator essentially 
obeys the same sort of kernel bounds in the periodic setting as it does in 
the non-periodic setting, and so the argument goes through as before. 
We omit the details. 



Remark 8.7. Theorem 8.2 asserts, roughly speaking, that if the en- 



ergy of the data is small in a large ball, then the energy will remain 



small in a slightly smaller ball for future times T; similarly, (67) as- 
serts that if the energy of the data is small outside a ball, then the 
energy will remain small outside a slightly larger ball for future times 
T. Unfortunately, this estimate is not of major use for the purposes of 



establishing Theorem 1.20, because energy is a supercritical quantity 
for the Navier-Stokes equation, and so smallness of energy (local or 
global) is not a particularly powerf ul co nclusion. To achieve this goal, 
we will need a variant of Theorem 8.2 in which the energy ^ / |mP is 
replaced by the enstrophy \ f |u;p, which is subcritical and thus able 
to control the regularity of solutions effectively. 
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Remark 8.8. It should be possible to extend Theorem 8^ to certain 
classes of weak solutions, such as mild solutions or Leray-Hopf solu- 
tions, perhaps after assuming some additional regularity on the solution 
u. We will not pursue these matters here. 



9. Bounded total speed 



Let (m, p, Mo, /, T) be an almost smooth finite energy solution. Applying 



the Leray projection P to dSl) (and using Corollary 4.3), we see that 



dtu = Au + PB{u, u) + Pf (68) 



for almost all times t, where B{u,v) = 0{'V{uv)) was defined in (12). 



As all expressions here are tempered distributions, we thus have the 



Duhamel formula (11), which we rewrite here as 



u{t) = e'^uo + / e(*-*')^(PC(V(MM)) + Pf){t') dt' . (69) 
Jo 



One can then insert the a priori bounds from Lemma 8.1 into m'M to 



obtain further a priori bounds on u in terms of the energy E{uq, /, T) 



(although, given that (53) was supercritical with respect to scaling, any 
further bounds obtained by this scheme must be similarly supercriti- 
cal). 

Many such bounds of this type already exist in the literature. For 
instance^ 



One can bound the vorticity a; := V x u in LfL\. norm [9], |36j : 
One can bound V^m in Lf^x'°° [2], [33] : 



4 4 

CK + l ' T CK + 1 ' 



More generally, for any a > 1, one can bound V"m in L^^ Lx 

m, IH]; 

For any A; > 0, one can bound t^d^u m L^^ [S]; 
One can bound Vu in L^J'^L^ [T6] ; 



For any r > and > 1, one can bound D[V*m in L^'"^'^'' ^ 

m, m-. 



2m 



For any 1 < m < oo, one can bound co in Lj"™ ^ L^"^ 
One can bound moments of wave-number like quantities |llj . 



^^These bounds are usually localised in both time and space, or are restricted 
to the periodic setting, and some bounds were only established in the model case 
/ = 0; some of these bounds also apply to weaker notions of solution than classical 
solutions. For the purposes of this exposition we will not detail these technicalities. 
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In this section we present another a priori bound which will be ab- 
solutely crucial for our localisation arguments, and which (somewhat 
surprisingly) does not appear to be previously in the literature: 

Proposition 9.1 (Bounded total speed). Let {u,p,uo, f,T) be a finite 
energy almost smooth solution. Then we have 

ll«llLjLsoao,T]xR3) <i?(«o,/,T)i/2Ti/^ + i?(«o,/,T). (70) 



We observe that the estimate (70) is dimensionally consistent with 



respect to the scaling (31). Indeed, if L denotes a length scale, then T 



scales like L^, u scales like L ^, and Eq scales like L, so both sides of 



(70) have the scaling of L. 



Before we prove this proposition rigorously, let us first analyse the 



equation (68) heuristically, using Littlewood-Paley projections, to get 



some feel of what kind of a priori estimates one can hope to establish 



purely from (68) and (53). For simplicity we shall assume / = 



for the sake of exposition. We consider a high-frequency component 



un '■= Pnu of the velocity field u for some N ^ 1. Applying P/v to (68), 
and using the ellipticity of A to adopt the heuristicj^ P/v^ ~ —N^Pn 
and P/vPV ~ NPn, we arrive at the heuristic equation 

dtUN = -N\n + 0{NPn{u^)). 

Let us cheat even further and pretend that Piy{u^) is analogous to u^un- 
(in practice, there will be more terms than this, but let us assume this 
oversimplification for the sake of discussion). Then we have 

dtUN = -N\n + 0{Nu%). 

Heuristically, this suggests that the high-frequency component un should 
quickly damp itself out into nothingness if ImatI ^ A^, but can exhibit 
nonlinear behaviour when \ui\i\ ^ A^. Thus, as a heuristic, one can 
pretend that un has magnitude ^ A^ on the regions where it is non- 
negligible. 



This heuristic, coupled with the energy bound (53 ), already can be used 



to informally justify many of the known a priori bounds on Navier- 



Stokes solutions. In particular, projecting (53) to the un component, 
one expects that 

IWnWl^li < (71) 

(dropping the dependencies of constants on parameters such as Eq, and 
being vague about the spacetime region on which the norms are being 
evaluated), which by Bernstein's inequality implies that 

^^One can informally justify this heuristic by inspecting the symbols of the 
Fourier multipliers appearing in these expressions. 
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However, with the heuristic that \un\ ^ on the support of un, we 
expect that 

summing in (and ignoring the logarithmic divergence that results, 
which can in principle be recovered by using Bessel's inequality to im- 
prove upon ([7l|)), we obtain a non-rigorous derivation of Proposition 

EH 



We now turn to the formal proof of Proposition |9.1[ All spacetime 

norms are understood to be over the region [0,T] x (and all spa- 
tial norms over R'^) unless otherwise indicated. We abbreviate Eq : = 
E{uo, f,T). From (53) and ^ we have the bounds 

hhrLi < (72) 

llVwIki^. < (73) 

\\M\Ll^mL\Li<ET- (74) 



We expand out u using (69). For the free term c^^Uq, one has by (18) 



\e'^noh^<t-'/'\\uo\\Li 



for t G [0,T], so this contribution to (70) is acceptable by (74). In a 
similar spirit, we have 

iie(*-*')Ap^(t')ii^^^ < {t-tr'/'\\pf{t')u. < {t-tr'^w)h. 

and so this contribution is also acceptable by the Minkowski and Young 



inequalities and (74). 



It remains to show that 



Jo 



By Littlewood-Paley decomposition, the triangle inequality, and Minkowski's 
inequality, we can bound the left-hand side by 



~E r /V^e(*-*')^O(PV(H(0)IUs° dt'dt. 
^ Jo Jo 

Using (27), and bounding the first order operator PV by A^ on the 
range of Pn, we may bound this by 

J2 [ [ exp{-c{t-t')N^)N\\PNO{uu){t')\\L^ dt'dt 

AT Jo Jo 



< 
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for some c > 0; interchanging integrals and evaluating the t integral, 
this becomes 

< ^ r N-'\\PNO{uu)it')\\L^ dt'. (75) 

We now apply the Littlewood-Paley trichotomy (see Section |2]) and 
symmetry to write 

PnO{uu) = X] PnO{un^un2) + PnO{un^un2) 

Ni^N N2<N Ni>N N2~Ni 

where u^- '■= Pnu. For Ni,N2 in the first sum, we use Bernstein's 
inequality to estimate 

\\PNO{U]\f-^U]\f^)\\Lcx, < lluArJI^oo ||UiV2 ||Lg° 

< iV(iV2/iVi)'/'l|Vn^JUil|Vn^JU|. 

For A^^i, N2 in the second sum, we use Bernstein's inequality in a slightly 
different way to estimate 

\\PNOiuN,ur,,)\\L^ < N^Oiur,,ur,,)\\Ll 
< N^\\unAl2J\un2\\li 
<N{N/m)'\\VuNA\Li\\Vur,,\\Li. 



Applying these bounds, we can estimate (75) by 

-T 



^ E E E (N^/my/' ! iivn^,(t')iiLiiiv^^,(t')iiLi df 

N Ni^N N2<N 

+ E E E w^i)' riiv«^.(oiUiiiv«^,(t')iu. dt'. 

N Ni>N N2'--Ni 

Performing the summation first, then using Cauchy-Schwarz, one 
can bound this by 

^ E E {N2/Niy/^aNiaN2+ E «^i"^2, 

Ni>lN2<Ni Ni>lN2r^Ni 

where 

ttN := ||Vm7v||l2l2- 



But from (73) and Bessel's inequality (or the Plancherel theorem) one 
has 



N 



and the claim (70) then follows from Schur's test (or Young's inequal- 

ity). 
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Remark 9.2. An inspection of the above argument reveals that the L 



oo 

X 

0,00 \ 



norm in (70) can be strengthened to a Besov norm defined by 

N 



Remark 9.3. An inspection of the proof of Proposition 9.1 reveals that 



the time-dependent factor T^^'^ on the right-hand side of Proposition 
9.1 was only necessary in order to bound the linear components 

e'^uo+ re(*-*')^(P/)(t') dt' 
Jo 



of the Duhamel formula (69). If one had some other means to bound 
these components in LlL"^ by a bound independent of T (for instance, 
if one had some further control on the decay of uq and /, such as L], 
and LlLl bounds), then this would lead to a similarly time- independent 



bound in Proposition |9]T| which could be useful for analysis of the long- 
time asymptotics of Navier-Stokes solutions (which is not our primary 
concern here). 

Remark 9.4. It is worth comparing the (supercritical) control given by 



Proposition 9.1 with the well-known (critical) Prodi-Serrin-Ladyzhenskaya 
regularity condition [351 SOI [271 CH 112] , a special case of which (roughly 
speaking) asserts that smooth solutions to the Navier-Stokes system 
can be continued as long as u is bounded in L'fL'^, and the equally 
well known (and also critical) regularity condition of Beale, Kato, and 
Majda |2j, which asserts that smooth solutions can be continued as 
long as the vorticity 

oj ■.= V xu (76) 

stays bounded in L\L°^ . 

Remark 9.5. As pointed out by the anonymous referee, one can also 
obtain L\L^ bounds on the velocity field m by a Gaghardo-Nirenberg 
type interpolation between the L/ bound on Vm from [16] with 



the LlL^^ bound on u arising from the energy inequality and Sobolev 
embedding. 



Although we will not need it in this paper. Proposition 9J^ when com- 
bined with the Picard well-posedness theorem for ODE yields the fol- 
lowing immediate corollary, which may be of use in future applications: 

Corollary 9.6 (Existence of material coordinates). Let {u,p,uq, f,T) 
be a finite energy smooth solution. Then there exists a unique smooth 
map $ : [0, T] x such that 

$(0, x) = X 

for all X G R^, and 

dt<^{t,x) = M($(t,x)) 
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zs a 



for all {t,x) G [0,r] x R^, and furthermore : — )■ R' 

diffeomorphism for all t G [0, T]. Finally, one has 

x)-x\< E{uo, /, T) + E{uo, /, T) 

for all (t,x) G [0,T] x R^. 

Remark 9.7. One can extend the results in this section to the periodic 
case, as long as one assumes normalised pressure and imposes the ad- 
ditional condition T < L^, which roughly speaking ensures that the 
periodic heat kernel behaves enough like its non-periodic counterpart 



that estimates such as ( 18 ) are maintained; we omit the details. (With- 



out normalised pressure, the Galilean invariance (33) shows that one 
cannot hope to bound the LjL'^ norm of u by the initial data, and even 
energy estimates do not work any more.) When the inequality T < 
fails, one can still obtain estimates (but with weaker bounds) by using 
the crude observation that a solution which is periodic with period L, 
is also periodic with period kL for any positive integer k, and choosing 
k to be the first integer so that T < (kL)'^. 

10. Enstrophy localisation 

The purpose of this section is to establish a subcritical analogue of 
Theorem |8.2[ in which the energy | J |Mp is replaced by the enstrophy 
I / Iwp. Because the latter quantity is not conserved, we will need a 
smallness condition on the initial local enstrophy; however, the initial 
global enstrophy is allowed to be arbitrarily large (or even infinite). 

Theorem 10.1 (Enstrophy localisation). Let {u,p,uo, f,T) be a finite 
energy almost smooth solution. Let B{xo, R) be a ball such that 

\\^o\\L2(Bixo,K)) + l|V X f\\LlLl{[0,T]xB{xo,R)) < ^ (77) 

for some 6 > 0, where := V x uq is the initial vorticity. Assume the 
smallness condition 

6^T + 6^E{uo,f,Ty/^T <c (78) 

for some sufficiently small absolute constant c > (independent of all 
parameters). Let < r < R/2 be a quantity such that 

r > C{E{uo, /, T) + E{uo, /, Ty/^T'/' + 6-') (79) 

for some sufficiently large absolute constant C (again independent of 
all parameters) . Then 

Remark 10.2. Once again, this theorem is dimensionally consistent 



(and so one could use (31) to normalise one of the non-dimensionless 



parameters above to equal 1 if desired). Indeed, if L is a unit of length. 
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then u has the units of L~^, u has the units of L~^, E{uo, f,T),r, R 
have the units of L, T has the units of L^, and 5 has the units of L~^/^ 
(so in particular 5^T and 6^E(uo, f,Ty^'^T are dimensionless) . 

Remark 10.3. The smallness of also comes up, not coincidentally, as 
a condition in the local well-posedness theory for the Navier-Stokes at 
the level of H^; see (46). The smallness of d^E{uo, f,Ty^'^T is a more 



artificial condition, and it is possible that a more careful argument 
would eliminate it, but we will not need to do so for our applications. 
For future reference, it will be important to note the fact that 6 is 
permitted to be large in the above theorem, so long as the time T is 
small. 

Remark 10.4. A variant to Theorem 110.11 can also be deduced from 
the resultp^ in [4,, Theorem D]. Here, instead of assuming a small 
condition on the enstrophy, one needs to assume smallness of quantities 
such as /j^3 dx for all sufficiently large Xq, and then regularity 

results are obtained outside of a sufficiently large ball in spacetime. 

We now prove the theorem. Let {u,p,uo, f,T), B{xo,R), 6,r be as 



in the theorem. We may use spatial translation symmetry (29) to 
normalise xq = 0. We assume c > is a sufficiently small absolute 
constant, and then assume C > is a sufficiently large constant (de- 
pending on c). We abbreviate Eq := E{uo, f,T). 

In principle, this is a subcritical problem, because the local enstro- 
phy ^ Jb{xo R) I'^l^ regularised versions thereof) is subcritical with 
respect to scaling ( [3l| ). As such, standard energy methods should in 
principle suffice to keep the enstrophy small for small times (using the 



smallness condition (78), of course). The main difficulty is that the lo- 
cal enstrophy is not fully coercive: it controls u (and, to a lesser extent, 
u) inside B{xo,R), but not outside B{xo,R); while we do have some 
global control of the solution thanks to the energy estimate (Lemma 
8.1[ ), this is supercritical and thus needs to be used sparingly. We 



will therefore expend a fair amount of effort to prevent our estimates 
from "leaking" outside B{xq,R)] in particular, one has to avoid the 
use of non-local singular integrals (such as the Leray projection or the 
Biot-Savart law) and work instead with more local techniques such as 
integration by parts. This will inevitably lead to some factors that blow 
up as one approaches the boundary of B{xo,R) (actually, for techni- 
cal reasons, we will be using a slightly smaller ball B{xo, R'{t)) as our 
domain). It turns out, however, that thanks to a moderate amount of 
harmonic analysis, these boundary factors can (barely) be controlled 
if one chooses exactly the right type of weight function to define the 
local enstrophy (it has to be Lipschitz continuous, but no better). 



19 



We thank the anonymous referee for this observation. 
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We turn to the details. We will need an auxiliary initial radius R' = 
i?'(0) in the interval [R — r/A, R] which we will choose later (by a pi- 
geonholing argument). Given this R', we then define a time-dependent 
radius function 

1 /■* 

R'{t) ■=R!-- ||n(s)|Uoo(R3) ds. 
c Jo 



From Proposition 9.1 one has 



R!{t) >R'- 0,{Eo + e;/'T^/^) 



and thus (by (79)) one has 

R'{t) >R-r/2 



if the constant C in (79 ) is sufficiently large depending on c. The reason 
we introduce this rapidly shrinking radius is that we intend to "outrun" 
all difficulties caused by the transport component of the Navier-Stokes 
equation when we deploy the energy method. Note that the bounded 



total speed property (Proposition 9.1) prevents us from running the 



radius down to zero when we do this. 

We introduce a time-varying Lipschitz continuous cutoff function 

ri{t, x) = min(max(0, c-^-^S^{R'{t) - \x\)), 1). 
This function is supported on the ball B{0,R'(t)) and equals one on 



B{0, _R'(t) — c°'^5~^), and is radially decreasing; in particular, from (79), 
we see that rj is supported on B{0, R) and equals 1 on B{0, R—r) if C is 
large enough. As t increases, this cutoff shrinks at speed ^||M(t)||ioo(R3), 
leading to the useful pointwise estimate 

dtr]{t,x) < -^\\u{t)\\L^^^3)\V,r]{t,x)\ (80) 

which we will use later in this argument to control transport-like terms 
in the energy estimate (or more precisely, the enstrophy estimate). 

Remark 10.5. It will be important that r] is Lipschitz continuous but no 
better; Lipschitz is the minimal regularity for which one can still control 
the heat flux term (see Is below), but is also the maximal regularity 
for which there is enough coercivity to control the nonlinear term (see 
Yq below). The argument is in fact remarkably delicate, necessitating a 
careful application of harmonic analysis techniques (and in particular, 
a Whitney decomposition of the ball). 



We introduce the localised enstrophy 

1 
2 



W{t) := - / \u{t,x)\'^r]{t,x) dx. 
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From the hypothesis (77) one has the initial condition 

W{0) < 6' (82) 

and to obtain the proposition, it will suffice to show that 

W{t) <, 6' (83) 

for all t E [0,T]. 



As u is almost smooth, W is C}. As in Section [sj we will compute the 
derivative dtW. We first take the curl of ^ to obtain the well-known 

vorticity equation 

dtu + {u ■ V)uj = Auj + 0{ujVu) + V x /. (84) 
This leads to the enstrophy equation 

dtllul"^ + {u ■ V)^|a;p = A^hul"^) - \Vuj\'^ + 0{ujujVu) + a; ■ (V x /). 

All terms in this equation are smooth. Integrating this equation against 
the Lipschitz, compactly supported 77 and integrating by parts as in 
Section |8] (interpreting derivatives of 77 in a distributional sense), we 
conclude that 

dtW = -Y, - F2 + I3 + n + n + n (85) 
where Yi is the dissipation term 



Y2 is the recession term 



2 

Y3 is the heat flux term 



Fi := / \Vu\'r], 

'R3 



Y, := -I I \u\'dtv 



R3 



2 

14 is the transport term 

1 
2 

Y5 is the forcing term 



R3 



F4 := - / \u\'^u ■ Vr], 

R3 



7r3 



and Yq is the nonlinear term 



Yq := / 0{ujLuVu)r]. 

JR3 
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The term Yi is non-negative, and will be needed to control some of the 



other terms. The term Y2 is also non-negative; by (80) we see that 

/ |u;nVr/| <c||n(t)|U^(R3)r2. (J 

We skip the heat flux term Y3 for now and use (86) to bound the 
transport term 1^ by 



;6) 



^7) 



Now we turn to the forcing term Y5. By Cauchy-Schwarz and (81) we 
have 

ini < w'/Mt) 



where 



Note from (77) that 



a{t) := ||V X f\\Li{B{o,R))- 



a{t) dt < 5. 



We return now to the heat flux term Y3. Computing the distributional 
Laplaciaij^ of t] in polar coordinates, we see that 

Y, < bit) 

where h{t) = bjii{t) is the quantity 

b{t) := c-°-^52^2 f 1^^^^ R'(t)a)\^ da+c-^'-^d^ [ \uj{t, x)\Ux 

Js^ J R'{t)-cO-'^5-^<\x\<R'{t) 

and da is surface measure on the unit sphere S"^. (Note that while 
At) also has a component on the sphere |a;| = R'{t) — c°'^(5~^, this 
component is negative and thus can be discarded.) 

To control b{t), we take advantage of the freedom to choose R'. From 
Fubini's theorem and a change of variables, we see that 



R 

R-r/4 Jo 



bn'it) dtdW < c"°-^52 f [ 

Jo Jb 



\co{t, x)\'^ dx. 



R3 

2JP /^0.1^ 



From Lemma 8.1, the right-hand side is 0(6 Eq/c ). Thus, by the 
pigeonhole principle, we may select a radius R' such that 

b{t) dt < 



^'^Alternatively, if one wishes to avoid distributions, one can regularise by a 
small epsilon parameter to become smooth, compute the Laplacian of the regu- 
larised term, and take limits as epsilon goes to zero. One can also rescale either R 
or S (but not both) to equal 1 to simplify the computations. 
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and in particular by ( 79 ) 

r-T 



/ b{t)dt<5' (89) 

^0 



if C is large enough. 



Henceforth we fix R' so that (89) holds. We now turn to the most 
difficult term, namely the nonlinear term Yq. Morally speaking, the 
Vu term in Yq has the "same strength" as u, and so Ig is heuristically 
as strong as 

A standard Whitney decomposition of the support of t], followed by 
rescaled versions of the Sobolev inequality, bounds this latter expres- 
sion by 



oaf M'vY'^j |Vu;pr/)). 

JR3 JR3 



If we could similarly bound Yq by this expression by an analogous argu- 
ment, this would greatly simplify the argument below. Unfortunately, 
the relationship between Vm and u is rather delicate (especially when 
working relative to the weight r/), and we have to perform a much more 
involved analysis (though still ultimately one which is inspired by the 
above argument). 

We turn to the details. We fix t and work in the domain 

n ■= B{0,R'{t)). 

We apply a Whitney-type decomposition, covering by a boundedly 
overlapping collection of balls Bi = B{xi,ri) with radius 

ri := ^mm{dist{xi,dQ),c°-^/6^). 
In particular, we have 

T] ~ c-^-^5\ (90) 
on B[xi, lOvi). We can then bound 

\Ye\<c-'-'6'y2n [ \co\'\Vu\. 

The first step is to convert Vm into an expression that only involves u 
(modulo lower order terms), while staying inside the domain Vt. To do 
this, we first observe from the divergence-free nature of u that 

An = VxVxu = Vxa;. 

Let ipi be a smooth cutoff to the ball S-Bj := -B(xj,3rj) that equals 1 
on 2Bi := B{xi,2ri). On 2Bi, we thus have the local Biot-Savart law 
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where v is harmonic on 2Bi. In particular, from Sobolev embedding 
one has 

||^^IUi(2i?.) < ll^i^^lL«/5(R,3) + \\u\\lI{2B,)- 

From Holder's inequality one has 

while from the mean value principle for harmonic functions one has 

\\Vv\\L'^^B,)<rf'^\\v\\L2i2B,y 

We conclude that 

\Wv\\L^(B,)<rf'^\\uj\\Ll(2B^)+rf'^\\u\\Lli2B,) 

and we thus have the pointwise estimate 

|Vm| < \V/\'^V{i)^UJ)\ + r~^''^\\i0\\Ll{2B,) + U^'^\\u\\lI(2B,) 

on Bi. We can thus bound l^el < ^6 i + ^6 2, where 



(91) 



and 
and 



:= |VA-V(^iu;)| +r, ^'^\\u\\lI(^2B,) 



Y.,2<C~-H'Y.-f"\\uU,2B, [ 

Let us first deal with 1^6.2, which is the only term that is not locally 
controlled by the vorticity alone. If the ball Bi is contained in the 
annnular region 

{xeVL:\x\> R'{t) 

which is the region where t] is not constant, then we use Holder to 
bound 



r 



-3/ 



2 I 



u\\lI{2B,) < ||m||l-(R3) 



and observe that = {Vrjl on B^. Thus, by (86), the contribution 

of this term to 16,2 is 0{c^'^Y2). If instead the ball Bi intersects the 
ball B{ Q, R' {t) — c^-^6~'^), then ~ c°-^5"^ and r/ ~ 1 on Bi, and we use 
Lemma 8.1 to bound '^^^||M||L2(2Bi) < c'^-^^S^eI^"^ , and then by (81), 

(g the contribution of this case is 0{c-^-^^6^Ei^^W) = 0{c°-^^W/T), 
thus 

Y,,,<c'-'Y2 + c'-''W/T. 
Now we turn to Ig,!- From Plancherel's theorem we have 

and thus 

WFiWlliB,) < ||w||l|{2B,)- 
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From Holder's inequality we thus have 

\^ ^ -0.1 c-2 3/2,, ||2 II II 

>6,i<c b 2^r/ ||a;||^6(B,)||w||L|{2B,)- 

i 

To deal with this, we let Wi denote the averages 
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Wi 



1 



135, 



1/2 



then 



^ 3/2 

Also, from the Sobolev inequality one has 



< II V(u;V^i) 11^2(^3) 

< ||Vcj||l2(3B,) + '"rlwIUKSB,) 

^ IIV7 II , 1/2 

< l|Vc^||L2(3i?,) 



and thus 



n,i < c-°-^5^5:rfz..||Vc.||i.(3^^) + c-°-^5^5^r,V. (92) 



To deal with the first term of (|92j), observe from (81) and (90) that 



(93) 



and in particular 



< c'-^^6-'W'/^rr^ (94) 
for all i. We may thus bound 



which by (90) and the bounded overlap of the Bi is 

Jn 

The second term of (92), c^'^'^S'^J^i^f'^iy is trickier to handle. Call 
a ball "large" if its radius is at least 10~^c~°'^(5~^ (say), and "small" 
otherwise. To deal with the small balls we use the Poincare inequality. 
From this inequality, we see in particular that 



135, 



\UJ 



2x1/2 



i\ J3B 



\3B 



ltd 



2N1/2 



j\ J3B 



< {rr^ / iV^n^^ 

JlOBi 
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whenever Bj intersect. (Indeed, the Poincare inequahty imphes 
that both terms in the left-hand side are within 0{{r~^ J^^^ |Va;p)"^/^) 
of l^iis^ fioB '^1') other words, we have 

\wr-W3\<r-"\l iVc^n^/^ (95) 



whenever Bi,Bj intersect. 



Now for any small ball Bi, we may assign a "parent" ball -Bp(i) which 
touches the ball but has radius at least 1.001 (say) as large as that of 
Bi. We may iterate this until we reach a large ball Ba(i), and write 

Wi < Wa(i) + ^ \Wpk(^i) - Wpk+^i)\ 
k>0 

where the sum is over all k for which p'^+^(i) is well-defined; note that 
this inequality also holds for large balls if we set a{i) = i. Taking cubes 
and using Holder's inequality, we obtain 



..3 ^ „..3 , Y^/i , L\10|„,, „,, |3 



fe>0 

and so we can bound c^^-^S"^ f^t'wf by 

< c-°-^5^ 5^r>^(, + c~'-'6^ 5^(1 + krY.'t\ 

i k>0 i 

If one fixes a large ball Bj, one easily checks that Ylira{i)=j ~ ^j^ 
thus 

i jr:rj>10-4cOl(5-2 

applying ( [94| ) and ( [93| ) we thus have 

i j 

Similarly, if one fixes a small ball Bj, one verifies that 

fe>0 i:pk(i)=j 

and thus 

fc>0 j j:rj<10-*c0i5-2 



From (94) (once) and (95) (twice) one has 
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and so we may bound the preceding expression by 
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< 



JlOB, 



which by (90) and the bounded overlap of the Bj can be bounded by 

Jn 

Putting the Yq^i bounds together, we conclude that 
collecting the bounds for Yi, . . . , we thus have 

To solve this differential inequality we use the continuity method. Sup- 
pose that < T' < T is a time for which 

sup W{t) < c-° °^5l (96) 

ie[0,T'] 

Then, if c is small enough, we can absorb the 0(c° °^5~^W^^/^Yi) term 
by the — Yi term, and can also use this bound and (78) to obtain 

-0. 



C -1553^^3/2 < ^-0.155^4^^ < c0.75y^/y 



and 

We thus have 



a{t)W 



1/2 < 



-0.005 



5a{t). 



dtW < c'^-^^W/T + c-°-°°^5a(t) + b{t). 



From Gronwall's inequality and (82), (88), (89) we thus have 



sup W{t) < c-°-°°^52_ 

ie[0,T'] 

For c a small enough absolute constant, this is (slightly) better than 
the hypothesis (96), and so from the continuity method (and (82)) we 
conclude that 



sup W{t) 



< 



-0.005 c-2 



te[o,T] 



and the claim (83 ) follows. The proof of Theorem 10.1 is now complete. 

Remark 10.6. As with Remark |8.4 we may adapt the proof of Theo- 
rem 10.1 to an annulus, replacing the ball B{xq,R) with an annulus 



B{xo, R')\B{xo, R) for some < R < R' with < r < R/2, {R'-R)/2, 
and replacing the smaller ball B{xq,R — r) with the smaller annulus 
B{xo, R' — r)\B{xo, R + r). To do this, one has to replace the cutoff 
T) (which was shrinking inside the ball B{xq, R) towards B{xq, R — r)) 
with a slightly more complicated cutoff (which is shrinking inside the 
annulus B{xo, R')\B{xo, R) towards the smaller annulus B{xo,R' — 
r)\B{xo, R + r)). However, aside from this detail, the proof method is 
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essentially identical and is omitted. Sending R' to infinity and using 
the monotone convergence theorem, we may in fact replace the annu- 
lus B{xo, R')\B{xo, R) with the exterior region Il^\B{xo, R), and the 
annulus B{xq, R! — r)\B{xo, R + r) with R.^\B{xq, R + r). 



Theorem 10.1 asserts, roughly speaking, that if the norm of the 
data is small on a ball, then for a quantitative amount of later time, 
the Hi norm of the solution remains small on a slightly smaller ball. 
As the norm is subcritical, we expect this sort of result to persist 
to higher regularities, in the spirit of [39]. It is therefore unsurprising 
that this is indeed the case: 

Proposition 10.7 (Higher regularity). Let {u,p,uo, f,T) be a finite 
energy almost smooth solution with T < T^,. Let B{xq,R), rj, 6, r 



obey the conditions (77), (78), ([79j) from Theorem 10.1. Then for any 
compact subset K in the interior of B{xo,R — r) and any k > 1, one 
can bound 

W'^'^uWl^ Ll{[0,T]xK) + l|V''"^^u||^2^2([o_r]xA') ^k,K,E{uoJ,T),5,n,R,Ak 1 

where 

k 

^k '■= ^ ||V-'Mo||L2(B(xo,i?)) + \\^^ f\\LfLl{[Q,T]xB{xo,R))■ 
j=0 

In particular, one has 

\\'^\\x'^([0,T]xK) ^k,K,E{uoJ,T),S,n,R,Ak 1- 

Proof. We allow all implied constants to depend on k, K, E{uo, /, T), 6, T^, R, A^. 
We introduce a compact set 

K C KiC K2C K3C K4C K5C B{xo, R-r) 

which each set lying in the interior if the next set. Let be a smooth 
function supported on K2 that equals 1 on f^i; we allow implied con- 
stants to depend on rj. 



We begin with the k = 1 case. From Theorem |10.1| one already has 

\\^\\LfLl{[0,T]xKr) + || Vw || ^2^2 ([o_r] xi^i) ^ 1- 

To pass from u to u, we use integration by parts. Since w = V x u and 
u is divergence-free, a standard integration by parts shows that 

\j \uj\^ri= [ \Vu\^7]+ f 0{\u\^V^7]). 
2 Jrs Jrs Jrs 

By Lemma 8.1, the error term is 0(1), and so we have 

|Vnp<l. 

K 
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Similarly, by replacing uj and u by their derivatives we also see that 

2 Jr3 Jr3 Jll3 



By Lemma 8.1, the error term is 0(1) after integration in time, and so 
we also have 

/ |V^m|^ dxdt < 1 
Jo Jk 

as desired. 



We now turn to the k = 2 case. This is the most difficult, as we cur- 
rently only control regularities that are half a derivative better than 

1 /2 

the critical regularity (which would place u in Hx ), and wish to boost 
this to three halves of a derivative above critical; this requires at least 
two iterations of the Duhamel formula. The arguments will be anal- 



ogous to the regularity arguments in Theorem 5.1 or Lemma 5.5 By 



(68) we see that urj obeys the truncated equation 

dt{r]u) - A{riu) = r]0{PV{uu)) + r]Pf + C(VuVr/) + 0{uV'^r]) (97) 

for almost all t. Meanwhile, from the k = 1 case and Lemma |8.1 we 
already have the estimates 

lkllLj°°L2([0,T]xR3) + \\^u\\L^L2^([Q^T\xKi) + II V^-u||2,2i2([o_T]xX4) ~ ^ (^8) 

and from the definition of A2 we have 

l|V%o||L2(S(xo,i?)) + \\^^ f\\L^Ll{[0,T]xB{xo,R)) H 1 (99) 

for J = 0,1,2. 



We claim that all terms on the right hand side of (97) have an L^L^([0, T] x 

norm of 0(1). The only difficult term here is r]PO{W{uu))\ the 
other three terms on the right-hand side are easily estimated in L^Ll. 



(and even in L^L^) using (98) and (99). We now estimate 



\\riO{PViuu) 



Il^L2([o,t]xR3)- 



We split uu = fjuu + {l — vi)uu, where is a smooth cutoff supported on 
K^^ that equals 1 on K^. For the contribution of the nonlocal portion 
{l — fj), one can use the smoothness of the kernel of the operator P away 
from the origin to bound this contribution by < ||C'('WM)||i4^i([o,T]xR3)) 



which is acceptable by (98); for future reference we note that this 
argument bounds this contribution on L^Ll. norm as well as in L^Ll. 
norm. For the local portion fjuu^ we discard the rj and P projections 
and bound this by 



LfL2([0,T]xR3)- 



But this is acceptable by (24). 
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We have now placed the right-hand side of (97) in LfLl{[0,T] x R'^) 
with norm 0(1). Meanwhile, from (99) the initial data UoT] is in ifJ(R^) 
with norm 0(1). Applying the energy estimate (23), we conclude that 



\ur]\ 



LfH: 



3/2- 



'([0,T]xR3) 



for any cr > 0. A similar argument (shifting the compact sets) also 
gives 



\urj 



3/2- 



+ \\U7] 



'([0,T]xR3) ' II""/ \\lIhI/'^-''{[Q,T]xH?) 

where rj' is a smooth function supported on that equals 1 on K^. In 
particular, by Sobolev embedding, on [0,T] x i^4, u is in L'^Ll^, Vm is 

in L'^Ll^nL^Lf^^, and V^m is in L^lI^^^, which together with ([98]) and 



< 1 



the Holder inequality now allows one to conclude that 0{W{'fiuu)) has 
an L'^H^{[0,T] x R^) of 0(1). Repeating the previous arguments, we 
now conclude that the right-hand side of (97) lies in L^Hl{[0,T] x R'^) 
with norm 0(1), and hence by (22) 



\\V'^\\l^H^{10,T]x-R3) + \\V'^\\l^^H3([0,T]x-R^) 

which gives the k = 2 case. 



The higher k cases are proven by similar arguments, but are easier 
as we now have enough regularity to place u in L^L!^([0,T] x K^) 
with norm 0(1); we leave the details to the reader. (For instance, to 
establish the k = 3 case, one can verify using the estimates already 
obtained from the k = 2 case that the right-hand side of (97) has an 
L2ffi([0, T] X R3) norm of 0(1). □ 



Remark 10.8. As in Remark 9.7 one can extend the results here to the 
periodic setting so long as one has T < and R < L; we omit the 
details. 



For our application to constructing Leray-Hopf weak solutions, we will 



need a generalisation of Theorem 10.1 to the case when one has hy- 



perdissipation. More precisely, we introduce a small hyperdissipation 
parameter e > 0, and consider solutions {u^^\p^^\uo, /, T) to the regu- 
larised Navier-Stokes equation, which are defined precisely as with the 
usual concept of a Navier-Stokes solution, but with ^ replaced by the 
regularised variant 

dtu'^'^ + (u(") • V)u(") = An(") - eA^M^") - Vp^'^ + f. (100) 

With hyperdissipation, the global regularity problem becomes much 
easier (the energy is now subcritical rather than supercritical), and 
indeed it is not difficult to use energy methods (see e.g. [32j) to show 
the existence of a unique almost smooth finite energy solution to this 
regularised equation {u^^\p^^\uo, f,T) from any given smooth finite 
energy data {uq, /, T). The energy estimate in Lemma 8.1 remains true 
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in this case (uniformly in e), and one easily verifies that one obtains 
an additional estimate 

e[ I \V'^u{t,x)\^ dtdx<E{uQj,T) (101) 
Jo Jn? 

in this hyperdissipative setting. One can also verify (with a some te- 



dious effort) that Proposition 9.1 also holds in this hyperdissipative 



setting as long as e is sufficiently small, basically because the hyper- 
dissipative heat operators e*^^"^'^ ^ obey essentially the same estimates 
(18), (27) as e*^ if < t < T and e is sufficiently small depending on 



T; we omit the details. 

One can define the vorticity a;'-^-* := V x u^^^ of a regularised solution 
as before. This vorticity obeys an equation almost identical to (84), 
but with an additional hyperdissipative term — eV^w''^^ on the right- 
hand side. One can then repeat the proof of Theorem 10.1 with this 
additional term. Integrating by parts a large number of times, one ob- 



tains a similar decomposition to (85) for the derivative of the localised 



e /j^3 |V^a;p?7 on 



enstrophy, but with the addition of a negative term 
the right-hand side, plus some boundary terms which are bounded by 
6(t), where 



E 

r=R'{t),R'{t)- 



ec-^-'5'R' / \Vu{t,ra)\' da 



52 



-F£c-°-'(5M \Vu{t,x)\^ dx 

J R'{t)-cO-'^5-^<\x\<R'{t) 

is a hyperdissipative analogue of b{t). By using the same averaging 
argument used to bound b{t) dt for typical R', one can also simulta- 
neously obtain a comparable bound for J^^ b{t) dt (taking advantage of 
the additional estimate (101)). The rest of the argument in Theorem 



10.1 works with essentially no changes; we omit the details. The proof 



of Proposition 10.7 is also essentially identical, after one notes that 



energy estimates such as (22) continue to hold in the hyperdissipative 



setting. Summarising, we we obtain 



Proposition 10.9. Theorem 10.1 and Proposition 10.7 continue to 



hold in the presence of hyperdissipation, uniformly in the limit e — > 0. 



11. Consequences of enstrophy localisation 



We now give a number of applications of the enstrophy localisation 
result. Theorem 10.1 Many of these applications resemble existing 



results in the literature, but with weaker decay hypotheses on the initial 
data and solution (in particular, we will usually only assumes either 
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finite energy or finite norm); the main point is that the locahsation 



afforded by Theorem 10.1 can significantly reduce the need to assume 



any stronger decay hypotheses. 

We begin with the observation that finite energy smooth solutions au- 
tomatically have bounded enstrophy if the initial data has bounded 
enstrophy: 

Corollary 11.1 (Bounded enstrophy). Let {u,p,uo, f,T) be an almost 
smooth, finite energy solution, such that the initial data (uq,/, T) has 
finite norm. Then u G ^^([0, T] x R^); in particular, {u,p, Uq, /, T) 
is an solution. 



Proof. Let 5 > be small enough (depending on E{uq, f,T),T) that 



the condition (78) holds. As {uQ,f,T) has finite norm, we have 

||Wo||l2(R3) + ||V X /||l1l2([0,T]xR3) < oo. 

By the monotone convergence theorem, we thus have for R sufficiently 
large that 

||l^o||L2{R3\B(0,i?)) + ||V X /||LiL2([o,T]x(R3\B{0,i?))) < ^• 



Applying Theorem 10.1 (inverted as in Remark 10.6), we conclude that 

\\^\\L^L'i{[0,T]x(R^\B(0,R+r))) + II Vw|| ^,2^2 qq.T] x {R3\_B(0,_R+r))) ~ ^ 

for some finite radius r, if R is sufficiently large; in particular, u lies in 
L'^Ll n L^H^ in the exterior region [0, T] x (R^\B{0, R + r)). On the 
other hand, as u is almost smooth, u also lies in L'^L'^ fl L^H^ in the 
interior region [0, T] x B{0, R + r + 1) (say). Gluing these two bounds 
together, we conclude that 

ueLrLln^Hli[0,T]xR'); 



meanwhile, from Lemma [8.1 [ one has 

ueLrLln^Hli[o,T]xR'). 

Since u is divergence-free and u = V x u, the claim then follows from 
Fourier analysis. □ 



Remark 11.2. From Corollary 5.8 we know that smooth solutions to 



the Navier-Stokes solutions can be continued in time as long as the 



norm remains bounded. However, Corollary 11.1 certainly does 
not allow one to solve the global regularity problem for Navier-Stokes, 
because the proof heavily relies on the solution u being complete rather 
than incomplete, thus it is (almost) smooth all the way up to the final 



time T, and not just smooth on [0,T). Instead, what Corollary 11.1 
does is to show that the solution from data is well-behaved when 
one is sufficiently close to spatial infinity; in particular, it does not 
prevent turbulent behaviour in bounded regions of spacetime. 
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Remark 11.3. If {u,p, uo, 0, T) is an almost smooth homogeneous finite 



energy solution, then by Lemma 8.1 we see that u(t) G if^(R^) for 



almost every time t G [0,T]. Applying the time translation symmetry 



(30) for a small time shift t^, we can then convert the finite energy 
data to data, and then by Corollary |ll.l| we conclude that in fact 
u{t) G Hl(R,^) for all non-zero times t G (0,T], and furthermore that 
u{t) is bounded in as soon as t is bounded away from zero. 



Since almost smooth solutions with normalised pressure are au- 
tomatically mild solutions, for which uniqueness was established 



in Theorem 5.4, we thus have uniqueness in the almost smooth finite 



energy category from smooth data: 

Corollary 11.4 (Unconditional uniqueness). Let (uq,/, T) be smooth 
data. Then there is at most one almost smooth finite energy solution 
{u,p,uo, f,T) with this data and with normalised pressure. 



This result resembles the standard "weak-strong uniqueness" results in 
the literature, such as those in |35], |10], [19], [20]. The main novelty 



here is the lack of decay hypotheses beyond the finite energy hypoth- 
esis; note that the almost smoothness of the solution gives plenty of 
integrability on compact regions of space, but does not imply any global 
integrability in space. 

Remark 11.5. We conjecture that one still retains uniqueness even if the 
data {uo, /, T^,) is merely smooth and finite energy, rather than smooth 



and H^. Note from Lemma 8.1 that u{t) has finite if^(R^) norm for 
almost every time t, which in principle allows one to enforce uniqueness 
after any given positive time (in the homogeneous case / = 0, at least), 
but it is not clear to the author how to prevent instantaneous failure 
of uniqueness at the initial time t = with only a smooth finite energy 
hypothesis on the initial data. It may however be possible to adapt the 
"weak-strong" uniqueness results of Germain [T^l |2D] to this category, 
perhaps in combination with the local control given by Theorem 

mm 



We now use the enstrophy localisation result to study solutions as they 
approach a (potential) blowup time T*. 

Proposition 11.6 (Uniform smoothness outside a ball). Let {u,p,uo, f,T~) 
be an incomplete almost smooth solution with normalised pressure 
for all times < T < T^. Then there exists a ball B{0, R) such that 

u,p,f,dtu G L-C,^([0,T,) X K) (102) 

for all k > and all compact subsets K of 'R.^\B{0, R) . 
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We remark that similar results were obtained in |1] assuming additional 
spatial decay hypotheses on the data at infinity, and in particular that 
/r3 l^lO, dx < oo. The main novelty in this proposition is that 

one only assumes square-integrability of Uq and its first derivatives, 
without any further decay assumption. 



Proof. From the argument in the proof of Corollary 11. 1| (noting that 
the bounds are uniform for all times T in a compact set), one can 
already find a ball -6(0, Rq) for which 

ueX\[^,%) X (R3\5(0,i?o))). 



Using Proposition 10.7 we then conclude the existence of a larger ball 
B{0,R) such that 

ueX\[0,T,) xK) 

for all A; > 1 and all compact subsets K of Ii^\B{0, R). From this, 
Sobolev embedding, and ^ (using the smoothness of the kernel of 



V'^A"^ away from the origin) we obtain (102) for u,p,f as desired. If 
one then applies ^ and solves for dtu one obtains the bound for dtu 
also. □ 



Remark 11.7. From (102) one can continuously extend u up to the 
portion {T*} x (R^\i?(0, R)) of the boundary (cf. the partial regularity 
theory in |lj). However, we were unable to demonstrate that u could be 
extended smoothly up to the boundary (or even that dtU is continuous 
in time at the boundary) . The problem is due to the non-local effects of 
pressure; the solution u could be blowing up at time T^, in the interior 
of B{0,R), leading (via ([9])) to time oscillations of the pressure in K 
(which cannot be directly damped out by the smoothness of the 
kernel, which only attenuates spatial oscillations) which by (|3| could 
lead to time oscillations of the solution u in K. Indeed, as Proposition 



1.12 shows, these time oscillations can have a non-trivial effect on the 



regularity of the solution. 



Remark 11.8. For future reference, we observe that Proposition [11.6 
did not require the full spacetime smoothness on /; it would suffice to 
have / G L^C^([0,T^,) x K) for all A; > and compact K in order 



to obtain the conclusion (102). This is because at no stage in the 



argument was it necessary to differentiate / in time. 



In a similar spirit, we may construct Leray-Hopf weak solutions that are 
spatially smooth outside of a ball for any fixed time T. More precisely, 
define a Leray-Hopf weak solution {u, p, uq, f, T) to smooth finite energy 
data {uq, /, T) to be a distributional solution u G X°([0, T] x R^) to ([3]) 
(after expressing this equation in divergence form) which is continuous 
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in time in the weak topology of L^(R'^), and which obeys the energy 
inequahty 

^ll«Wllii(R3) + l|Vn(t)||i.(i,,) dx < E{uo, f,T). (103) 

The existence of such solutions was famously demonstrated by Leray 
[31] for arbitrary finite energy data (mq, f,T); the singularities of these 
solutions were analysed in a vast number of papers, which are too 
numerous to cite here, but we will point out in particular the seminal 
work of Caffarelli, Kohn, and Nirenberg [3]. 

Our main regularity result for Leray-Hopf solutions is as follows. 

Proposition 11.9 (Existence of partially smooth Leray-Hopf weak 
solutions). Let {uq, f, T) he smooth data. Then there exists a Leray- 
Hopf weak solution {u,p,uo, f,T) to the given data and a ball B{0,R) 
such that u is spatially smooth in [0,T] x (R^\i?(0, i?)) (i.e. for each 
t e [0,T], u{t) is smooth outside of B{0, R)). 

Again, similar results were obtained in [1] under stronger decay hy- 
potheses on the initial data. We also remark that weak solutions which 
were only locally of finite energy, from data of uniformly locally finite 
energy, were constructed in [2S]; the ability to localise the weak solu- 
tion construction in this fashion is similar in spirit to the results in the 
above proposition. 



Proof. (Sketch) We use a standard hyperdissipatioEp^ regularisation 
argument. Let e > be a small parameter, and consider the al- 
most smooth finite-energy solution {u^'^\p^^\uo, f, T) to the regularised 
Navier-Stokes system (100), which can be shown to exist by energy 



methods. By Proposition 10.9 we can extend Theorem 10.1 and Propo- 



sition 10.7 (and thence Proposition 11.6), to these regularised solutions 
u'^^\ with bounds that are uniform in e as e — > 0. As a consequence, we 
can find a ball 5(0, R) independent of e such that for every compact set 
K outside of S(0, R) and every A; > 0, V'^m^") lies in L~L^([0, %] x K) 
uniformly in A^. If we then extract a weak limit point u of the u^'^\ 
then one by standard arguments one verifies that m is a Leray-Hopf 
weak solution which is spatially smooth outside of 5(0, R). □ 



21t 



It may also be possible to use other regularisation methods here, such as veloc- 
ity regularisation, to construct the Leray-Hopf weak solution; however, due to the 



delicate nature of the proof of the localised enstrophy estimate (Theorem lO.l I, we 
were not able to verify that this estimate remained true in the velocity-regularised 
setting, uniformly in the regularisation parameter, due to the less favourable vor- 
ticity equation in this setting. 
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Remark 11.10. As before, we are unable to demonstrate regularity of 
u in time due to potential non-local effects caused by the pressure, 
which could in principle cause singularities inside B{0,R) to create 
time singularities outside of B{0,R). 

Remark 11.11. Uniqueness of Leray-Hopf solutions remains a major 
unsolved problem, for which we have nothing new to contribute; in 
particular, we do not assert that all Leray-Hopf solutions from smooth 



data obey the conclusions of Proposition 11.9 However, if (mq, /, oo) 
is globally defined smooth data, the above argument gives a single 
global Leray-Hopf weak solution {u,p, Uq, f, oo) with the property that, 
for each finite time T < oo, there exists a radius Rt < oo such that u is 
smooth in [0,T] x (R^\_B(0, i?)). If we restrict to the case / = 0, then 



from (103) we see that || VM(t)||i2(R3) must become arbitrarily small 
along some sequence of times t = tn going to infinity. If ||Vu(t)||L2(j^3) 
is small enough depending on E{uo, 0, oo), then standard perturbation 
theory arguments (see e.g. [21]) allow one to obtain a smooth, bounded 
enstrophy solution from the data u{t) on {t, +oo), which by the unique- 
ness theory of Serrin |l0] must match the Leray-Hopf weak solution u 
on {t, -|-oo). As such, we conclude in the homogeneous smooth case 
that one can construct a global Leray-Hopf weak solution which is spa- 
tially smooth outside of a compact subset of spacetime [0, +oo) x R^. 
Again, we emphasise that this global weak solution need not be unique. 



12. Smooth solutions 



The purpose of this section is to establish Theorem 1.20[ iii). To do 
this, we will need the ability to localise smooth divergence-free vector 
fields, as follows. 

Lemma 12.1 (Localisation of divergence- free vector fields). LetT > 0, 
0<Ri<R2<R3< Ra, and let u : [0, T) x (5(0, Ra)\B{0, Ri)) 
be spatially smooth and divergence-free, and such that 

u,dtu e Lra^([0,T) X {B{0,Ra)\B{0,R,))) 

for all k >0 and 

/ u{t,x) -n da{x) = (104) 

J \x\=r 

for all Ri < r < Ri and t G [0,T), where n is the outward normal 
and da is surface measure. Then there exists a spatially smooth and 
divergence- free vector field u : [0,T) x {B{0, Ra)\B{0, Ri)) — > R'^ which 
agrees with u on [0,T) x (5(0, -R2)\-B(0, -Ri)), hut vanishes on [0,T) x 
(5(0, i?4)\i?(0, i?3)). Furthermore, we have 

u,dtueLrC^A[^,T) X (5(0,i?4)\5(0,i?i))) 

for all k>0. 
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Finally, if we have 

1 < 2i?2 < i?3 < R2 

then we have the more quantitative bound 

\\u\\l^H'={[0,T)x{B{0,R4)\B{(},Ri))) \\u\\l^h''+i([o^t)x(B(0,R4)\B{0,Ri))) 

(105) 

for any k. (This latter property will come in handy in the next section.) 



Note that the hypothesis (104) is necessary, as can be seen from Stokes' 



theorem. Lemmas of this type first appear in the work of Bogovskii [3]. 

Proof. One could obtain this lemma as a consequence of the machin- 
ery of compactly supported divergence-free wavelets |28], but for the 
convenience of the reader we give a self-contained proof here. 

Let X denote the vector space of all divergence-free smooth functions 
u : -8(0, i?4)\-B(0, Ri) — t- obeying the mean zero condition 

/ u{x) ■ n da{x) = (106) 

J \x\=r 

for all Ri < r < _R4, and such that \\u\\ck(^(o,R4)\B(o,Ri)) < C)0 for all k. 
It will suffice to construct a linear transformation P : X ^ X that is 
bounded^ from C''+^ to C^, i.e. 



l|-P^llc"={(0,i?4)\B(0,_Ri)) ^Ri,R2,R3,R4,k ll^llc'=+2((0,R4)\B(0,_Ri)) 

for all k > 0, and such that Pu equals u on B{0, R2)\B{0, Ri) and 
vanishes on _B(0, -R4)\-B(0, i^s), as one can then simply define u{t) : = 
Pu{t) for each t G [0,T). 

We now construct P. We work in polar coordinates x = ra with 
Ri ^ r ^ Ra and a G S*^ (thus avoiding the coordinate singularity at 
the origin), and decompose u{r,a) as the sum of a radial vector field 
Ur{r,a)a for some scalar field Ur, and an angular vector field Ua{r,a) 
which is orthogonal to a; thus, for fixed r, Ua{r) can be viewed as a 
smooth vector field on the unit sphere S"^ (i.e. a smooth section of 
the tangent bundle of S"^). The divergence-free condition on u in these 
coordinates then reads 

drUr{r) + -V« ■ uJr) = (107) 
r 



while the mean zero condition ( 106 ) reads 



Ur{r, a) da = 0. 

52 

^^One can reduce this loss of regularity by working in more robust spaces than 
the classical C'^ spaces, such as Sobolev spaces i?** or Holder spaces C'^'", but we 
will not need to do so here. 
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Note that either of these conditions imphes that drUr{r) has mean zero 
on S"^ for each r. From (107) and Hodge theory we see that 



rA^^V adrUrir) + vir) 



where inverts the Laplace-Beltrami operator A^^ on smooth mean 
zero functions on S"^, and v{r) is a smooth divergence-free vector field 
on S'^ that varies smoothly with r. 

Let Tj : [Ri^R^ -t- R+ be a smooth function that equals 1 on [Ri^Rq] 
and vanishes on [R^^R^. We define 

Ur := rj{r)ur 

and 

Ua{r) = rA'^VadrUrir) + r]{r)v{r) 

and 

Tu := u := UrQ + Ur,- 



One then easily verifies that u is smooth, divergence-free, obeys (106), 
depends linearly on u, equals u on B{0, R2)\B{0, Ri), and vanishes 
on B{0, Ri)\B{0, R4). It is also not difficult (using the fundamental 
solution of A~^) to see that T maps (7'^+^ to (with some room to 
spare). The claim follows. 



Finally, we prove (105). It suffices to show that 

\\Tu\\Hk(B{0,Ra)\B{0,R2)) 1 

whenever k > 0, and m G X is such that 

\\u\\h''+2(B(0,R4)\B{0,Ri)) ^ 1- 



Henceforth all spatial norms will be on -6(0, R3)\B{0, R2), and all im- 
plied constants may depend on k. As u has an H^~^^ norm of 0(1), Ur 
and hence Ur has an H'''^^ norm of 0(1) also. As for Ua, we observe 
from the Leibniz rule that 

Ua = VUa + {rdrri{r))A~^V aUr{r). 

As u has an H^^^ norm of 0(1), we have r~'^'W^^dlua has an norm 
of 0(1) whenever i + j < k + 1, which (using elliptic regularity in 
the angular variable) implies that r~*V^9^'UQ, has an norm of 0(1) 
whenever i + j < k. This gives m = + an norm of 0(1) as 
claimed. □ 



We can now establish Theorem 1.20[ iii): 



Theorem 12.2. Suppose Conjecture \l.S\ is true. Then Conjecture 1.1 S\ 
is true. 
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Proof. In view of Corollary 5.8 , it suffices to show that if (m, p, uq, f, ) 
is an incomplete mild solution up to time T^:, with uq, f spatially 



smooth in the sense of Conjecture 1.19 then u does not blow up in 
enstrophy norm, thus 

limsup ||M(t)||j/i(R3) < oo. 



Let -R > be a sufficiently large radius. By arguing as in Corollary 



11.1, we have 



ueLrHl{R\B{0,R)) 
and thus the blowup must be localised in space: 

limsup ||M(t)||^^i(B(o,i?)) < oo. 



(108) 



By Proposition 11.6 and Remark 11.8 (and increasing R if necessary) 
we also have 



u,pj,dtu G LrC^'([0,T.) X (5(0,5i?)\i?(0,2i?))) 



(109) 



for all > 0. From Stokes' theorem and the divergence-free nature of 
M, we also have 



u{t, x) ■ n da{x) 







|x|=r 



for all r > and t G [0, T). Applying Lemma [l2.1[ we can then find a 
spatially smooth divergence-free vector field u : [0, T) x (5(0, 5-R)\-B(0, 2R)) — )■ 

which agrees with u on B{0, 3R)\B{0, 2R) and vanishes outside of 
5(0, 4i?), with 

u, dtu G L,°^C^(5(0, 5i?)\5(0, 2R)) (110) 

for all A; > 0. We then extend u by zero outside of 5(0, 5/2) and by u 
inside of 5(0, 2i2), then u is now smooth on all of [0, T) x R^. 

Let ?7 be a smooth function supported on 5(0, 5i?) that equals 1 on 
5(0, 4i?). We define a new forcing term / : [0,T) x R^ — )■ R by the 
formula 

f := dtu + {u-\/)u- /\u + \/{pri), (111) 
then / is spatially smooth, supported on 5(0,5/2) and agrees with / 



on 5(0,3/?). From this and (110), (109) we easily verify that 

/GLr5i([0,r,)xR3). 



Note from taking divergences in (111) and using the compact support 
of prj, u, f that 



pr] 
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Thus, {u,pri, u{0), f, T~) is an incomplete pressure-normalised (and 
hence mild) solution with all components supported in -8(0, 5-R). If we 
then choose a period L larger than lOR, then we may embed -8(0, 5-R) 
inside R^/LZ^ and obtain an incomplete periodic smooth solution 
{l{u), L^prj), l{u{0)), L{f),T~ , L), where we use t(/) to denote the ex- 
tension by zero of a function / supported in -8(0, 5-R), after embedding 
the latter in H^/LZ^. By construction we then have 

i{f)eLrHl{[0,n)xRyLZ'). 

As {T*} has measure zero, we may arbitrarily extend / to [0, T ^,] x 
R^/-LZ^ while staying in Lf'H^. Applying either Conjecture 1.9 (and 
the uniqueness component to Theorem 5.1) or Conjecture 1.10, we 
conclude that 

l{u) e L'^Hl{[0,n) X RVLZ=^) 
which implies (since u and u agree on -8(0, R)) that 
u e L^Hl{[0,T,) X -8(0, i?)) 



which contradicts (108). The claim follows. 



□ 



Observe that if we omit the embedding of -8(0, 5-R) in R^/-LZ'^ in the 
above argument, we can also deduce Conjecture |1.19|from Conjecture 



1.18[ Since Conjecture |1. 19 clearly implies Conjecture |1. 18| as a special 
case, we obtain Theorem 1.20[ iii). 



Remark 12.3. The referee has pointed out a variant of the above argu- 
ment using the partial regularity theory of Caffarelli, Kohn, and Niren- 
berg j3], which allows one to partially reverse the above implications, 
and in particular deduce Conjecture [1.8 from Conjecture 1.191 We 



sketch the argument as follows. Assume Conjecture 1.19 and assume 



for contradiction that Conjecture |1.8| fails, thus there is a periodic solu- 
tion with smooth inhomogeneous data which first develops singularities 
at some finite time T, and in particular at some location (T, xq). We 
may extend the solution beyond this time as a weak solution. Applying 
a periodic version of the theory in |3] , we see that the set of singularities 
has zero one- dimensional parabolic measure, which among other things 
implies that the set of radii r > such that the solution is singular 
at (T, x) for some x with |x — Xo| = r has measure zero. Because of 
this, one can find radii r2 > ri > such that the solution is smooth 
in the annular region {{t,x) '■ < t < T;ri < \x — Xq\ < r2}. By 
smoothly truncating the solution u to this annulus as in the proof of 



Theorem 12.2 one can then create a non-periodic mild solution to 
the inhomogeneous Navier-Stokes equation with spatially smooth data 
which develops a singularity at (T, xq) while remaining smooth up to 



time T, contradicting Conjecture 1.19 (when combined with standard 



uniqueness and regularity results, such as those in Theorem 5.4) 
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13. Smooth finite energy solutions 



In this section we establish Theorem 1.20 v). It is trivial that Conjec- 



ture 1.14 implies Conjecture 1.13, so it suffices to establish 



Theorem 13.1. Suppose that Conjecture 1.13 is true. Then Conjee 



ture 1.14 is true. 



Proof. Let (mo,0,T) be smooth homogeneous finite energy data. Our 
task is to obtain an almost smooth finite energy solution {u,p, Uq, 0, T) 
with this data. We allow all implied constants to depend on uq. 

We use a regularisation argument. Let Nn be a sequence of frequencies 
going to infinity, and set u^^^ := Pkn^^Uq, then Mq"^ converges to Uq 
strongly in L^(R'^), and (mq"'', 0, T) is smooth data for each n. Thus, 
by hypothesis, we may find a sequence of almost smooth finite energy 
solutions , 0, T) with this data. 

One could try invoking weak compactness right now to extract a solu- 
tion, but as is well known, one only obtains a Leray-Hopf weak solution 
by doing so, which need not be smooth. So we will first work to es- 
tablish some additional regularity on the sequence (after passing to a 
subsequence as necessary) before extracting a weakly convergent limit. 

Since the (mq"''*,0,T) are uniformly bounded in energy, we see from 
Lemma 18.11 that 

lk("1xO([0,T]xR3) < 1. (112) 

Now let < To < T/2 be a small time. From ( 112[ ) and the pigeonhole 



principle, we may find a sequence of times t^^'G [0, tq] such that 

lk("Hr("))l|Hi(R3)<ro-^ 

Passing to a subsequence, we may assume that r*^") converges to a limit 
r G [0,ro]. If we then take r' G [t, 2ro] sufficiently close to r, we may 



apply Lemma 5.5 and conclude that 



(say) for all sufficiently large n. Passing to a further subsequence, we 
may then assume that u^'^\t') converges weakly in if^°(R^) (and thus 
locally strongly in H^) to a limit Mq G if^°(R'^). By hypothesis, we 
may thus find an almost smooth solution {u',p', u'q, 0,T — r') with 
this data. 



Meanwhile, by time translation symmetry (30), (u^'^\- + r'),p^^\- + 
t') , u^"\t') , 0,T — t') is also a sequence of almost smooth solutions. 
Since ^("'^(r') converges locally strongly in if^(R^) to u'q, we would 
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like to conclude that u^'^\t + r') also converges locally strongly to u{t) 
in if^(R'^), uniformly in t G [0,T — r']. This does not quite follow 
from the standard local well-posedness theory in Theorem 5^, because 
this theory requires strong convergence in the global if^(R^) norm. 
However, we may take advantage of the local enstrophy estimates to 
spatially localise the local well-posedness theory, as follows. 

Let e > be a small quantity (depending on the solution u' = {u',p', Uq, 0, T- 
r')) to be chosen later, let i? > be a sufficiently large radius (depend- 
ing on e and {u' ,p' ,Uq,0,T — r')) to be chosen later. Since u'q is in 
H^^(R^), we see from monotone convergence that 

||'«o||/fiO(R3\B(0,_R)) ^ ^ (113) 

if R is sufficiently large depending on e. Since the u*^"^(r') converge 
locally strongly in H^(Il^) to u'q, we conclude that 

\\u^'^\t')\\hio(^b{o,wr)\b{o,r)) ^ £ 



if n is sufficiently large depending on R, e. Applying Theorem |10.1 
conclude (if R is large enough depending u'q and T — r') that 



we 



||^(")(^. _^ r')||xi([0,T-r']x(B(0,9i?)\_B(0,2R))) ^ £ 

for n sufficiently large depending on R, e. Using Duhamel's formula 



(and Corollary 4.3) repeatedly as in the proof of Proposition 10.7 



we 



may in fact conclude that 

||9jM^"^(- + r') II 2,^00 ^6([o,T-r']x(B(O,8iJ)\B{O,3i?0)) ^-u'.T ^ (114) 

(say) for i = 0, 1, taking R large enough depending on u',T,e to en- 
sure that the contributions to the Duhamel formula coming outside 
B[0,9R) or inside B{0,2R) are negligible, and taking n sufficiently 
large as always. 



4.3 



We let p^^^ be the normalised pressure, defined by ([9]); by Corollary 
p^'^\t) and p^'^\t) differ by a constant C{t) for almost every t. Using 
(|9|, (114) and Lemma 8.1, we see that 



lb^"^l|Lj°°H|([0,T-T']x(B(0,7i?)\B(0,4il))) Hu' ,T 

if R is large enough depending on u',T,e. 



Applying Lemma 12.1 , we may find divergence-free smooth vector fields 



u 



[r', T] X R3 ^ R3 which agree with on [r', T] x 5(0, 5R) but 
vanish outside of [r',T] x B{0,6R), with 

||9j-U^"^(- + r')||2,J«H5{[0,r-r']x(B(0,8i?)\B(0,3R))) ^u' ,T £ (115) 

(say) for n sufficiently large and i = 0, 1. 
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Let rjhe a smooth function that equals 1 on B{0, 6R) and is supported 
on B{0,7R), and obeys the usual derivative bounds in between. We 
then consider the smooth solutions 

+ r'), r^p^-H- + r ), w^-n^'), /^"\ T - r') (116) 

where 

/(") := + . Vu^"^ - Au(") + V(r/p(")))(- + r') 

By construction, /' and /*^"^ are smooth and supported on [0, T — r'] 



X 



{B{0,7R)\B{0,5R)), and the (116) are smooth, compactly supported 
solutions. From the preceding bounds on u^"'\p^"-^ we see that 

||/^"^||L(°°//l{[0,T-r]xR3) <n',T £ 

for n sufficiently large. 



Also, using (113), (115) we have 



for n sufficiently large. If e is sufficiently small, we conclude from the 



local well-posedness theory (Theorem 5.4) that 



(• + — w'||xi([0,T-t']xR3) ^«',T £ 



and in particular 

+ r') - M'|Ui([0,T-r']xB(0,fl)) ^u',T ^ 

for n large enough. Sending e to zero (and R to infinity) we conclude 
that M^"^(- + r') converges weakly to u' . In particular, we see that any 
weak limit of the u'-"-' is smooth on [r', T] x (and furthermore, the 
weak limit is unique in this spacetime region). 

The above analysis was for a single choice of r. Choosing r to be a 
sequence of times going to zero (and repeatedly taking subsequences 
of the ■u^"^-' and diagonalising as necessary) we may thus arrive at a 
subsequence m*^"-* with the property that there is a unique weak limit u 
of the u*^"-*, which is smooth on (0, T] x R^. If we then set p by ([9]), we 
see on taking distributional limits that (m,p, uo,0,T) is a Leray-Hopf 
weak solution to the initial data (uo,0,T). 

To finish the argument, we need to show that (u,p, uo,0,T) is almost 
smooth at (0,a;o) for every xq G R^. Fix xq, and let > be a 
large radius. As Uq is smooth, ||%||//1(b(xo,5R)) is finite, and hence 

(re- 



m{B{xo,bR)) is uniformly bounded. Applying Theorem 



10.1 



calling that the -u*^") have uniformly bounded energy), we conclude (for 
R large enough) that there exists < r < T such that ||-u(") ||xi([o,r]xB(xo,4-R)) 
is uniformly bounded in n. Using Duhamel's formula as in Proposition 
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and noting that u^"^ is uniformly smooth on i?(xo,4i?), we con- 
clude that ||M''"^||L^c'=((o,r]xB(xo,3-R)) uniformly bounded for all k >0. 
Taking weak limits, we conclude that 

u G L^C'^((0,r] X B{xo,3R)) 



for all k >0. From this and ^ (and Lemma 8.1), we also see that 



peLrC''{{0,T]xB{xo,2R)) 

for all /c > 0. Using ([s]), we conclude that 

dtu e L~C^X(0,r] X B{xo,2R)) 

for all A; > 0. A similar argument also shows that 

e L^C^{{0,r] X B{xo,2R)) 

uniformly in n. From this, we see that the V^m'-"-' are uniformly Lips- 
chitz in a neighbourhood of (0,xo). Since V^m*-"^ converges weakly to 
the smooth function V^u in (0, T] x R^, and also converges strongly at 
time zero in HKR^) to the smooth function V^Mq, we conclude that 
V^M can be extended in a locally Lipschitz continuous manner from 
(0, T] X to [0, T] X R^ in such a way that it agrees with V^Mq at 
time zero. 

Now we consider derivatives V*^p of the pressure near (0, Xq). Let e > 
be arbitrary. Then by the monotone convergence theorem, we see that 
if i?' > is a sufficiently large radius, then 



and thus 

for n large enough. 



||^io|Ui(R3\B(xo,i?')) ^ ^ 

I (")ll <- 
\Uo \\l2{rs\B{xo,R')) ~ ^ 



By Theorem 8.2, we conclude that if R' is large enough, there exists a 
time < r < T such that 

||^i^"''lU?°L2([0,T]x(R3\B(xo,2iJ'))) ~ ^ 

and hence on taking weak limits 

||w||Lj°°L2([0,r]x{R3\B(xo,2R'))) ~ ^• 

On the other hand, as V^u is continuous at t = 0, u{t) converges in 
C^{B{xo,2R')) to Mo as t for any A; > 0. From this and (|9| (and 
the decay of derivatives of the kernel of A^^ away from the origin) we 
see that 

limsup \V^p{t,x) - VVo(a;o)| <fc e 

(t,x)-5-(0,xo);i>0 

for any > 0, where po is defined from uq using (|9|. Sending £ — > and 
i?' — i- oo we conclude that VV extends continuously to V^Pq{xq) at 
(0, xo), and thus extends continuously to V^o on all of the initial slice 
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{0} X R^. By ([3]) we conclude that dtV^u also extends continuously to 
the initial slice, with the Navier-Stokes equation ([s]) being obeyed both 
for times t > and times t = 0. We have thus constructed an almost 
smooth finite energy solution (-u,/), -Uq; 0, T) as desired. □ 



Remark 13.2. We emphasise that Theorem |13.1| only establishes exis 
tence of a smooth finite energy solution 
and not uniqueness; see Remark 11.5 



assuming Conjecture 1.13), 
However, it is not difficult to 
see from the argument that one can at least ensure that the solution 
constructed is independent of the choice of time T, and can thus be 
extended to a single global smooth finite energy solution. (Alterna- 
tively, from Lemma [8.1| we see that the enstrophy of the solution will 



become arbitrarily small for a sequence of times going to infinity, so 
for a sufficiently large time one can in fact construct a global smooth 
solution by standard perturbation theory techniques.) 

Remark 13.3. One can modify the above argument to also establish 



Conjecture 1.14 with a non-zero Schwartz forcing term /, provided of 



course that one also assumes Conjecture 1.13 can be extended to the 
same class of /. We have not however investigated the weakest class of 
forcing terms / for which the argument works, though certainly finite 
energy seems insufficient. 



14. Quantitative bounds 



In this section we prove Theorem 1.20[ vi). We begin with some easy 
implications. Firstly, it is trivial that Conjecture 1.17 implies Con- 



jecture 1.16, and from the local well-posedness and regularity theory 



Proposition 5.6) 



in Theorem 5.4| (or Corollary 5.8) we see that Conjecture 1.16 



plies Conjecture 1.15, which in turn implies Conjecture 1.13 (thanks to 



im- 



Next, we observe from Theorem 5.1 and Lemma 5.5 that given any 
data (mo,0,T), there exists a time < r < T such that one has an 
mild solution (u,p, uo,0,r) with u{t) smooth. If Conjecture 1.13 



holds, then one can then continue the solution in an almost smooth 
finite energy manner (and hence in an almost smooth manner, 
thanks to Corollary 11.1) from r up to T. Normalising the pressure 



of this latter solution using Lemma 4J^ and gluing the two solutions 
together, we obtain an mild solution up to time T. From this we 



see that Conjecture 1.13 implies Conjecture 1.15 



Now we show that Conjecture 1.16 implies Conjecture 1.17 
that one has homogeneous data (mo,0,T) with 



Suppose 



F0||/fi(R3) 



< A<oo. 
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By Conjecture 1.16 (which imphes Conjecture 1.15) we may obtain a 
mild solution {u,p,uo,0,T), which is smooth for positive times. 
Our objective is to show that 

lkllLj°°//i([0,T]xR3) 1- 



Let £ > be a quantity depending on A to be chosen later. We may 
assume that T is sufficiently large depending on e, A, otherwise the 



claim will follow immediately from Conjecture 1.16 Using Lemma 8.1 
and the pigeonhole principle, we may then find a time < Ti < T with 
Ti <A 1 such that 

l|VM(ri)|U2(R3) <£. 

Meanwhile, from energy estimates one has 

\\u{Ti)\\L2^n^) <A 1. 

On [Ti, T], we split u = ui + v, where ui is the linear solution ui{t) : = 
e^^~'^^\(Ti) and v := u ~ Ui. From (21) one thus has 



\Ul\\xo 



<A 1 



and 



|Vmi||xo < e. 



From (11), ( 22 ) one has 



||^^IUi([Ti,T]xR3) < \\0{uiVui +uiVv + vVui +fVt;)||i2i2([r^_r]xR3). 

We now estimate various contributions to the right-hand side. We 
begin with the nonlinear term 0{v'Wv). By Holder (and dropping the 
domain [Ti, T] x for brevity) followed by Lemma 8.1 we have 



^ X t X X X 



1/2 



1 1/2 



1 1/2 



^ II l|3/2|| II 
<- II l|3/2 



A similar argument gives 



^ ^II^IUi 



and 



and 



\0{uiVv)\\l2li < l|Vf |Uo||Vmi||5/o||mi||^/^^ 



<Ae'''\\v\\x^ 



LOCALISATION AND COMPACTNESS FOR NAVIER-STOKES 



83 



and thus 



I^IU^ <Ae'" + e'/'\\v\U. + \\v\\f,. 



If e is small enough depending on A, a continuity argument in the T 
variable then gives 

and thus 

ll'"l|xi{[Ti,T]) 1- 

Using this and the triangle inequality, we conclude that Conjecture 



1.16 implies Conjecture 1.17 



We now turn to the most difficult implication: 

Proposition 14.1 (Concentration compactness). Suppose that Con- 



jecture 1.15 is true. Then Conjecture 1.16 is true 



We now prove this proposition. The methods are essentially those of 
|17j (which are in turn based in [1], [IBJ), which treated the (more 
difficult) critical analogue of this implication; indeed, one can view 



Proposition 14.1 as a subcritical analogue of the critical result |T71 
Corollary 1]. For the convenience of the reader, though, we give a 
self-contained proof here, which does not need the full power of the 
machinery in the previously cited papers because we are now working 
in a subcritical regularity rather than a critical regularity such as 
if^/^, and as such one does not need to consider the role of the scaling 



symmetry (31 ) 



We first make the remark that to prove Conjecture |1.16[ it suffices to 
do so with the condition 

||wo||^fi(R3) < A (117) 

replaced by (say) 

||«o||hioo(r3) < A (118) 
To see this, observe that if we take data uq in if^(R^), then from 



Theorem |5.4| and Lemma [5~5] there exists a time Ti > depending only 
on A such that 

ll^l|Lj°°/i'i([0,min(T,Ti)]xR3) 1, 

and such that 

||m(Ti)||j/100(r,3) <A 1 



if T > Ti. From this and time translation sym metry (30 ) we see that we 
can deduce the if^(R^) version of Conjecture 1.16 from the H^^'^CR^) 
version. 

Now suppose for contradiction that the Hl'^^(R,^) version of Conjecture 



1.16 failed. Carefully negating the quantifiers, we can find a sequence 
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, 0, T*-")) of smooth homogeneous solutions, with T*^"^ 
uniformly bounded, and Mq"^ uniformly bounded in if^°°(R^), such that 

lim ||w^"^||l-^i([o,t(")]xr3) = oo- (119) 



By Lemma we may assume that these solutions have normalised 
pressure. 

If we were working on a compact domain, such as R^/Z^, we could now 

(n) 

extract a subsequence of the Mq that converged strongly in a lower 
regularity space, such as if^^(R^/Z^). But our domain R^ is non- 
compact, and in particular has the action of a non-compact symmetry 
group, namely the translation group t^^u^x) := u{x — Xq). However, 
as is well known, we have a substitute for compactness in this setting, 
namely concentration compactness. Specifically: 

Proposition 14.2 (Profile decomposition). Let Mq"'' G if^°°(R^) be a 

sequence with 

lim sup ||4"^||^^ioo(R3) < A, 

n—^oo 

and let e > 0. Then, after passing to a subsequence, then there exists 
a decomposition 

J 

(n) V"^ , (n) 

J=l ' 

where \ J\ <a,£ 1, Wi^, . . . ,Wjfi G if^°°(R^), xj"^ G R^, and the re- 
mainder rl^^ obeys the estimates 

limsup ||rJ"^||Hioo(R3) < A 

and 

lim sup ||rJ"^||Ljo(R3) < e. (120) 

Furthermore, for any 1 < J < j' < J, one has 

\x^;^ -xPl^oo, (121) 

and for any 1 < j < J, the sequence T_^(„)rl(^^ converges weakly in 
i/i°°(R3) to zero. 

Finally, if the Uq""* are divergence-free, then the wj^ and r'^'' are also 
divergence- free. 



Proof. See e.g. [18]. We sketch the (standard) proof as follows. If 
ll'^o"'' IUg°(R.^) — ^ sufficiently large n then there is nothing to 

prove (just take J = and Tq"'' := Uq"^). Otherwise, after pass- 
ing to a subsequence, we can find a sequence G R^ such that 
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l^o"^(^i"^)l — (say)- The sequence r_ (n)ui^^ is then bounded in 
H]f^(R^) and bounded away from zero at the origin; by passing to a fur- 
ther subsequence, we may assume that it converges weakly in if^'^°(R^) 
to a hmit Wi, which then has an if^°°(R^) norm of '^A,e 1 ^.nd is asymp- 
totically orthogonal in the Hilbert space if^''°(R'^) to T_jn)U^Q'\ We can 
then decompose 

(n) , (n),l 

and from an application of the cosine rule in the Hilbert space if^™(R^) 
one can verify that 

limsup ||'"o"^'"^||hioo{r3) < A'^ — c 

n— >oo 

for some c > depending only on e, A. We can then iterate this 
procedure Oj^si^) times to obtain the desired decomposition. □ 



We apply this proposition with a value of e > depending on A, T to be 
chosen later. The Wj^ lie in H]f^(R.^), and thus by the assumption that 



Conjecture 1.15 is true, we can find mild solutions {wj,pj, wj^, 0, T) 



with this data. By Theorem 5.1 we have 



II 11x100 < OO 

for each 1 < j < J, and to abbreviate the notation we adopt the 
convention that the spacetime domain is understood to be [0, T] x R^. 



Next, we consider the remainder term rl^\ From (21) one has 

l|e*M"^IUioo<A 



while from (120) one has 



for n sufficiently large. Interpolating between the two, we soon conclude 
that 

l|e*M"^IUi <A,T 

for some absolute constant c > 0. If we take e sufficiently small de- 
pending on A, T, we can use stability of the zero solution (see Theorem 



5.1 one could also have used here the results from [6]) to conclude the 
existence of a mild solution {r^^\pi"\r^\o, T) with this data, with 
the estimates 

l|r(")||xi <A,T^^; (122) 
from Theorem 15.11 we then also have 

lk^"^IUioo <A,T 1. 

We now form the solution 
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where the velocity field ■u'-"^ is given by 

i=i ' 

the pressure field p*^"^ is given by ([9]), and the forcing term /*^") is given 
by the formula 

/(") := - Au^") - P5(u("\ 

This is clearly a mild solution, with 

11x100 <A,T,e 1- 

We now estimate f''^\ From (68) for the solutions r {n)Wj + r*^"'^ we 

^3 



have an expansion of /'•"^ purely involving nonlinear interaction terms: 

/W= ^ PO(V(r^,„)t.„r^(„)«;,0) 
i<j</<j ' ' 

i<i<J 

In particular, from the triangle inequality and translation invariance 
we have 

ll/^"^llL?Li < E ll^(VK-,r^(;)_,W«^,'))llL?Li 

i<i<i'<J 

+ ll^(VK,r_^wr(")))|U?L- 



i<i<J 



But by (121) and Sobolev embedding, r ,(n)_ and r_ xnjr*^"-' are 

bounded in L^L^ and converge locally uniformly to zero, and so we 
conclude that 



lim ||/(")|L2^. =0. 



From this and the stability theory in Theorem 5.4| we conclude that 



for n large enough, there is an mild solution 
with 

lim = 0, 

n— >oo 

and in particular 

limsup ||M^"^|U-jyi([o,T]xR3) < oo. 



By the uniqueness theory in Theorem 5.4, this solution must agree with 



the original solutions mJ^, T^")) on [0,T(")] x R^; but then 



we contradict (119). Proposition 14.1 follows. 
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15. NON-EXISTENCE OF SMOOTH SOLUTIONS 



In this section we establish Theorem 1.12 Informally, the reason for the 



irregularity is as follows. Assuming normalised pressure, one concludes 
from Q that 

If one then differentiates this twice in time, using ^ to convert time 
derivatives of u into Am plus lower order terms, integration by parts 
to redistribute derivatives, we eventually obtain (formally, at least) a 
formula of the form 

d^p = 0{A'^V\AuAu)) + l.o.t.. 

But if Uq is merely assumed to be smooth and in H^, then Au can grow 
arbitrarily fast at infinity at time t = 0, and this should cause p to fail 
to be at time zero. 

We turn to the details. To eliminate the normalised pressure assump- 
tion, we will work with Vp instead of p, thus we will seek to establish 
bad behaviour for Vdfp at time t = 0. For technical reasons it is con- 
venient to work in the weak topology in space. The key quantitative 
step is the following: 

Proposition 15.1 (Quantitative failure of regularity). Let Uq : — )■ 

be smooth, divergence-free, and compactly supported, and let ip : 
R^ — 7- R &e smooth, compactly supported, and have total mass ip = 
1. Let R,M,e > 0. Then there exists a smooth divergence-free com- 
pactly supported function ui which vanishes on B{0,R) with 

\\ui\\hi(r^) ~ ^ 

and such that if {u, p,uo + ui,0,T) is a mild ( and hence smooth, by 



Proposition \5.d^ solution with data {uq + Ui,0,T), then 

Vdfp{0, x)ijj{x) dx 



> M. (123) 



Let us assume this proposition for now and conclude Theorem 1.12 
We will use an argument reminiscent of that used to establish the 
Baire category theorem or the uniform boundedness principle. Let 
^ : R^ — )■ R be a fixed smooth, compactly supported function with 
total mass 1. We will need a rapidly decreasing sequence 

e«>£(2)>...>0 

of small quantities to be chosen later, with each e^'^^ sufficiently s mall 
depending on the previous e^^\ . . . Applying Proposition 15.1 



recursively starting with uq = 0, one can then find a sequence of 
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smooth, divergence- free, and compactly supported functions for 
n = 1, 2, . . . such that 



H IIhi(r3) ^ e'- 

with uS"^ vanishing on 5(0, l/e^")), such that if 0' ^^"^) 

is a mild (and hence smooth) solution with data 

(n) (1) , , (n) 

then 

V9jV"H0,a;)^/'(x) dx| > l/e^"^ (124) 

Furthermore, each depends only on e^^\ . . . , and in particular 
is independent of 

By the triangle inequality (and assuming the e^*^-* decay fast enough), 
the data is strongly convergent in H^(R^) to a limit Uq = "'''"^ 



i/i , with 



If we make each sufficienly small depending on ul^\ and hence 

on e^^\ . . . ,e^"\ then the m^"^ will have disjoint supports; as each -u^"^ 
is smooth and divergence- free, this implies that uq = ^i"'' is also 

smooth and divergence-free. 



Applying Theorem 5.1, we may then take the times T*^"^ = 1 (if the 
e^"^ are small enough), and {u''^\p^"'\u''^\o, 1) will converge to a mild 
solution (m,p, Mo,0, 1) in the sense that m^"-* converges strongly in 
X^([0, 1] X R^) to u. Indeed, from the Lipschitz stability property we 
see (if the e^"'^ decay fast enough) that 

Ik -w(")|Ui{[0,l]xR3) 



Also, are bounded in X^{[0,1] x R^) by 0(1). Using (g and 

Sobolev embedding, this implies that 



and so if one sets 



and 



b-p^"^IUr^i([o,i]xR3)<^^"+'\ 



F(")(t) := / Vp^''\t,x)ij{x) dx 



F{t) := / Vp{t,x)ilj{x) dx 

iR3 

then from integration by parts we have 



F-F^-'>U^no,,^.<e^-+'\ (125) 
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Meanwhile, each F^"^^ is smooth, and F continuous, from Proposition 



5.6, and from (124) one has 

|52^(n)(0)| > 1/5 W. 

In particular, if is sufficiently small depending on F^"^^ (which 

in turn depends on e^^\ . . . one has from Taylor's theorem with 

remainder that 

|ir(n)(2(e("+i))0.i) _ 2FW((e("+i))0-i) + F(")(0)| ^ („) 



Applying (125), we conclude that 

|^(2(£("+i))0-i) - 2F((£("+i))0-i) + F(0)| 



(^(n+l))0.2 

if £("+1) is sufficiently small depending on €^^\ . . . In particular, 

\F{2h)-2F{h) + F{0)\ 
hmsup ^-^ — — ^— = +00 

which by Taylor's theorem with remainder implies that F is not smooth 
at 0. 

We claim that the data mq gives the desired counterexample to Theo- 
1.12[ Indeed, suppose for contradiction that there was a smooth 



rem 



solution Mq, 0, T) for some T > 0. By shrinking T we may assume 



T < 1. By Lemma 4.1 we see that p{t) has normalised pressure up 
to a constant for almost every t, and thus after adjusting p{t) by that 
constant, Mq, 0, T) is a mild solution. Using the uniqueness 



property in Theorem 5.1, we conclude that u = u, and p{t) and p{t) 
differ by a constant for almost every t, and hence (by continuity of both 
p and p) for every t. In particular, Vp = Vp, and so 

F{t) = / Vp{t,x)ilj{x) dx. 

But as p is smooth on [0,T] x R^, F is smooth at 0, a contradiction. 

Remark 15.2. The above argument showed that Vp failed to be smooth 
at t = 0; by using ^ we conclude that the velocity field u must then 
also be non-smooth at t = (though the velocity u has one more degree 
of time regularity than the pressure p). Thus the failure of regularity 
is not just an artefact of pressure normalisation. Using the vorticity 



equation (84) one can then show a similar failure of time regularity 
for the vorticity, although again one gains an additional degree of time 
differentiability over the velocity u. 



The irregularities in time stem from the unbounded growth of high 
derivatives of the initial data. If one assumes that all spatial derivatives 
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of Mo are in L^(R^), i.e. that uq € H°°(R^) then one can prove itera- 
tiveljj^ that all time derivatives of u and p at time zero are bounded, 
and also have first spatial derivatives in H'^(R^) (basically because 
the first derivative of the kernel of the Leray projection is integrable at 
infinity). In particular, u and p now remain smooth at time 0. 



It remains to establish Proposition 15.1 Fix uo, ipy e, and let ui 



be a smooth divergence-free compactly supported function ui vanishing 
on B{0,R) with if^(R^) norm 0{e) to be chosen later. Let {u,p, up + 



Ml, 0, T) be a mild solution with this given data. By Theorem 5.1 
this is a smooth solution, with all derivatives of u,p lying in L'^L'i. 
From Lemma 14.11 we thus have 

Vp = -VA-^didj{uiUj) (126) 

for almost all times t. But both sides are smooth in [0,T] x R^, so 
this formula is valid for all times t (and in particular at t = 0). In 
particular, we may apply a Leray projection P to ([s]) and conclude 
that 

dtu = Au + PB{u,u). (127) 



We differentiate (126) once in time to obtain 

Vdtp = -2VA~^didj{uidtUj). 
Expanding out dtUj using (|3]), we obtain 

Vdtp = -2VA-^ didjiuiAuj) + 0{A-^V^{uPB{u,u))). 
Writing 

didjiuiAuj) = -2didjiidkU^)idkUj)) + OiV\uu)) 
we thus have 

Vdtp = 2VA-^didj{dkUidkUj)+0{A~^V^{uu))+0{A-^V%uPB{u,u))). 



We differentiate this in time again and use (127) to obtain 

Vd'^p = 4:VA-^didj{dkUidkAuj) 

+ 0{A-^V%(yu)VPB{u,u))) 
+ 0{A-^V''{udtu) 
+ 0{A-^V%{dtu)PB{u,u))) 
+ 0{A-^V%uPB{u,dtu))). 



'^^We thank Richard Melrose for this observation. 
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We can write dkUidkAuj = —{Aui){Auj) + C(V(VmAm)), so that 

V(9,V = -4:VA-^didj{AuAuj) 
+ 0{A-^V\VuAu) 
+ 0{A-^V%{Vu)VPB{u,u))) 
+ 0{A'^V\udtu) 
+ 0{A-'V\{dtu)PB{u,u))) 
+ 0{A-^V%uPB{u,dtu))). 
Integrating this against ip, we may thus expand 

5 

/ V9fp(0, x)ij{x) dx = 4Xo + V 0{Xi) 



where 



Xo := [ {didjV A-^i))AuiAuj 

Xi ■= [ (V^A-V)VwAm 

X2 := /" (V^A-V)VnVP5(u,u) 

Xs:= [ (V^A-V)«<9tM 

X^:= [ {V^A-^i)){dtu)PB{u,u) 

JRS 

X5:= ! {\/^A-'^lj)uPB{u,dtu), 
with all expressions being evaluated at time 0. 



From (127) and Sobolev embedding one has 

||5tM(0)||i2(R3) <„0 1+ lkl||/f2(R3). 

Meanwhile, if e is small enough, we see that 

||M(0)||/fi(R3) So 1 

and thus from the Gagliardo-Nirenberg inequality 

I|m(0)||l^(r3) So (1 + lkl||H2(R3))^/^ 

From many applications of the Sobolev and Holder inequalities (and, 
in the case of X5, an integration by parts to move the derivative off of 
dtu), we conclude that 

\Xi\ So,^ (1 + ll'"l||/f|(R3))^''^ 
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for i = 1,2, 3, 4, 5. In a similar spirit, one has 



+ I|wi||h|(r3)) 



To demonstrate (123), it thus suffices to exhibit a sequence : R'^ — )■ 
of smooth divergence-free compactly supported vector fields sup- 
ported outside of B{0, R) such that 



L,(")||2 



with 
and 



|mS"^I|hi(r3) 



I (")ll 



We construct u^^^ explicitly as the "wave packet" 

where 61,62,63 is the standard basis, Xq G R^ is a point (independent 
of n) outside of B{0,R + 1) to be chosen later, and 

'^^'^\x) = x{^) sin(n^ ■ x)rj 

where ^ G R^ is a non-zero frequency (independent of n) to be chosen 
later, 77 G R"^ is a non-zero direction, and x • R-^ — )■ R is a smooth 
bump function supported on -8(0, 1) to be chosen later. Note from 

construction that u['^^ is smooth, divergence-free, and supported on 
B{xo,l), and thus vanishing on B{0,R) for Rq > R + 1. One can 
compute that 



and 



("■)ll -1/2 
^1 lki(R3)<X^ 



as long as x is not identically zero. To conclude the theorem, it thus 
suffices to show that 



Jr3 

if Rq and n are large enough. 



n 



Observe that 

uf'\x) := n~^/^sin(n,^ ■ {x — Xq))x{x — Xq){^ x rj) + 0{n^^^'^) 
and similarly 

Au^^\x) := sin« ■ {x - Xo))x{x - Xo){^ x r/) + 0(^-^/2) 
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and so by choosing % appropriately, and using the Riemann-Lebesgue 
lemma, it suffices to find xo,C,V ^ such that 

{didjVA-',p){^ X tjUC X v)Axo) + 0. 

But as ■?/' has mean one, we see that V^A~^-?/'(a;o) is not identically zero 
for xq large enough, and the claim follows. 
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